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ABSTRACT 


A  mathematical  model  is  proposed  for  a  solid- liquid 
catalytic  reaction  carried  out  in  a  nonideal  continuous -flow- 
stirred-tank  reactor  (CFSTR).  The  model  considers  a  nonideal 
CFSTR  as  a  combination  of  a  perfectly  backmixed  reactor  and  a 
by-pass  region  and  describes  the  effect  of  external  mass  transfer , 
intraparticle  diffusion  and  chemical  reaction  (first  order,  irre¬ 
versible)  and  the  nonideality  in  the  reactor  operation  on  the 
overall  reaction  rate.  The  derived  rate  expression  permits  eval¬ 
uation  of  the  various  parameters  which  characterize  the  individual 
phenomena. 

Attempts  were  made  to  verify  the  model  experimentally. 

The  experimental  set  up  consisted  of  a  glass  CFSTR  and  the  related 
accessories.  The  reaction  employed  was  the  hydrolysis  of  ethyl 
acetate  catalysed  by  a  cation  exchange  resin.  The  experimental 
work  was  performed  in  three  parts. 

The  first  part  involved  the  determination  of  the  response 
of  the  reactor  to  a  step  change  in  the  feed.  This  was  done  to 
evaluate  the  extent  of  nonideality  at  different  agitator  speeds. 
The  results  indicate  that  the  reactor  exhibited  no  nonideality 
in  operation  under  the  experimental  conditions  employed. 


. 


The  second  part  was  concerned  with  the  evaluation,  of  the 
significance  of  the  external  film  diffusion  by  carrying  out  the 
reaction  at  various  agitation  rates.  The  problem  was,  however, 
complicated  due  to  settling  of  the  resin  particles  at  low  stirrer 
speeds.  For  all  stirrer  speeds  at  which  the  particles  were 
completely  in  the  bulk  phase,  the  external  film  mass  transfer 
resistance  was  negligible. 

In  the  last  and  more  extensive  phase  of  the  study  the  rate 
constant,  k,  and  the  effective  diffusivity,  D,  for  ethyl  acetate 
in  the  resin  phase  were  determined  at  different  temperatures  by 
carrying  out  the  reaction  at  high  stirrer  speeds.  Effectiveness 
factors  for  the  catalyst  were  calculated  for  different  temperatures 
and  particle  diameters.  The  results  indicate  that  the  chemical 
reaction  is  the  rate  limiting  step  for  the  hydrolysis  of  ethyl 
acetate  under  the  experimental  conditions  employed. 

These  results  are  discussed  in  light  of  the  work  done  by 
previous  investigators.  Though  only  a  partial  verification  of 
the  model  was  achieved,  significant  knowledge  was  gained  on 
various  aspects  of  the  system. 
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CHAPTER  I 


INTRODUCTION 

The  investigation  of  a  heterogeneous  catalytic  reaction 
is  complicated  by  the  fact  that  in  such  a  system  there  is  an  in¬ 
terplay  between  the  true  chemical  reaction  rate  and  the  physical 
transport  phenomenon.  Depending  on  the  reaction  system  and  the 
operating  conditions  either  one  may  contribute  the  dominating 
influence  on  the  overall  rate  process. 

It  is  postulated  that  in  the  presence  of  a  relative 
motion  between  a  solid  particle  and  the  ambient  fluid,  a  thin 
film  of  the  liquid  envelopes  the  particle.  The  reactants  diffuse 
through  this  film  and  then  into  the  catalytic  particle  where 
chemical  reaction  takes  place.  To  isolate  only  one  of  these 
rate  processes  for  study,  it  is  necessary  to  devise  conditions 
so  that  all  the  other  effects  play  no  role  in  the  observed 
kinetics.  However,  these  extremes  are  not  always  experimentally 
accessible.  Often  within  the  range  of  conditions  for  a  chemical 
reaction  the  individual  rate  processes  are  coupled  and  then  it 
is  necessary  to  account  for  their  combined  effects  on  the  reaction 


rate. 
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With  a  given  reaction  system,  the  effect  of  varying 
turbulence,  catalyst  size  and  temperature  in  an  experimental 
study  should  indicate,  respectively,  the  importance  of  external 
film  diffusion,  intra-particle  diffusion  and  chemical  reaction. 

It  is  convenient  to  study  these  aspects  in  a  Continuous-Flow- 
Stirred-  Tank  Reactor  (CFSTR).  The  experimental  data  do  not 
require  differentiation  or  integration,  nor  is  an  especially 
precise  analysis  required.  Unlike  tubular  reactors  the  intensity 
of  turbulence  is  independent  of  the  flow  rate  and  can  be  easily 
controlled  by  the  degree  of  agitation.  Therefore  the  effect  of 
film  mass  transfer  may  be  more  easily  evaluated  or  eliminated  in 
a  CFSTR.  Similarly,  the  temperature  gradients  within  the  bulk 
phase  can  be  eliminated. 

By  definition,  the  contents  of  an  ideal  CFSTR  are 
perfectly  mixed;  they  are  at  uniform  composition  and  temperature. 

At  steady  state  the  exit  stream  has  the  same  composition  as  the 
fluid  within  the  reactor.  However,  deviations  from  the  ideal 
behaviour  may  occur  at  low  stirrer  speeds  which  give  imperfect 
mixing.  This  nonideality  results  in  nonuniform  composition  in 
the  reactor  and  causes  changes  in  the  extent  of  the  reaction. 

The  particular  system  chosen  for  experimental  study  was 
the  hydrolysis  of  Ethyl  Acetate  using  a  suspension  of  acidic  ion 
exchange  resin  paricles  as  catalyst.  The  ester  hydrolysis  proceeds 
according  to  the  reaction: 
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C2H5COOCH3  +  H20  - —  CH3COOH  +  CgH^OH 

The  reaction  is  fairly  rapid  at  moderate  temperatures  in 
the  presence  of  the  catalyst  and  proceeds  at  a  negligible  rate 
without  the  catalyst.  By  employing  a  large  excess  of  water  the 
rate  of  the  reaction  can  be  approximated  by  the  linear  relation 

Rate  =  kC 

This,  in  addition,  ensures  low  concentrations  of  the 
products  in  the  resulting  mixture.  The  reaction  can  therefore 
be  assumed  as  virtually  irreversible. 

Previous  experiments  with  similar  systems  show  wide 
disagreement:  some  investigators  reported  significant  film- 

diffusion  effects  while  others  found  simultaneous  intra-particle- 
diffusion  and  chemical  reaction  as  the  rate  controlling  process. 
Again,  while  there  exists  considerable  literature  on  the  theoretical 
aspects  of  nonideal  behaviour  of  backmix  reactors,  some  workers  have 
reported  only  minor  deviations  even  in  the  absence  of  any  stirring. 

In  face  of  these  seemingly  divergent  observations,  the 
aim  of  this  investigation,  then,  was  to  establish  the  relative 
importance  of  different  rate  processes  under  varying  operating 
conditions  for  the  reaction  under  study,  to  ascertain  to  what 
extent  nonideality  in  flow  exists  and  also  to  investigate  a  means 
of  accounting  for  this  nonideality  in  prediction  of  the  actual 


conversion  in  the  reactor. 


CHAPTER  II 


THEORY 


II,-  A.  GENERAL 


A  catalyst  is  a  substance  which  by  virtue  of  its  presence 
in  the  reaction  mixture,  changes  the  rate  of  chemical  reactions 
but  is  not  consumed  stoichiometrically  in  the  reactions.  Catalysts 
offer  an  alternate  route  for  the  reactions  to  take  place.  This 
property  of  a  catalyst  is  advantageously  employed  in  reactions 
which  otherwise  require  extreme  operating  conditions. 

In  a  heterogeneous  system,  a  catalyst  is  instrumental  in 
the  chemical  conversion  in  the  sense  that  the  reaction  takes  place 
at  or  near  the  solid  surface.  Therefore  the  reactants  present  in 
the  ambient  fluid  have  to  be  transferred  to  the  interface  and  the 
products  away  from  it.  In  the  main  body  of  a  well-agitated  fluid, 
the  transport  takes  place  by  a  rapid  eddy  diffusion.  However,  due 
to  frictional  effects  near  the  solid  surface,  the  catalyst  particle 
is  surrounded  by  a  thin  film  of  the  liquid.  Mass  transfer  through 
this  relatively  stagnant  film  occurs  only  by  molecular  diffusion 
and  convection. 

It  is  recognized  that  in  a  heterogeneous  catalysic  reaction, 
the  following  steps  occur  successively: 
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(a)  The  reactants  diffuse  through  the  liquid  film  to 
the  particle  surface, 

(h)  The  reactants  then  diffuse  into  the  interior  of 
the  catalyst, 

(c)  Chemical  reaction  takes  place,  and 

(d)  The  products  diffuse  outwards  through  the 
surrounding  liquid  film  to  the  bulk  of  the  fluid. 

As  the  chemical  reaction  occurs  simultaneously  with  the 
intra-particle  diffusion,  they  provide  resistances  in  parallel 
and  they  both  offer  resistance  in  series  with  the  liquid  film  mass 
transfer  resistance.  Whenever  one  of  the  steps  offers  major 
resistance,  that  step  may  be  considered  as  the  rate  controlling 
step.  Resistances  of  the  different  steps  can  and  do  vary  greatly 
from  each  other.  Also  their  relative  importance  varies  with  the 
operating  conditions.  The  individual  resistances  and  their 
dependence  on  the  operating  conditions  must  be  known  in  order  to 
express  the  overall  rate  expression  in  terms  of  the  concentrations 
in  the  bulk  fluid  stream. 

II-B .  LIQJJID  FILM  DIFFUSION 

Many  mechanisms  have  been  postulated  to  explain  the  mass 
transfer  of  solute  from  the  bulk  fluid  to  the  particle  surface. 
That  this  transfer  of  material  is  very  closely  related  to  the 
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existing  flow  pattern  in  the  vicinity  of  the  particles  and  to  the 
fluid  properties  is  at  once  obvious.  The  present  qualitative 
picture  is  of  continuous  increase  in  turbulence  from  the  interface 
to  the  main  body  of  liquid.  The  concept  of  a  film  as  proposed  by 
Lewis  &  Whitman  (1924)  can  be  used  as  a  convenient  basis  for  a 
mental  picture  of  the  .localization  of  diffusional  resistance  in  a 
region  near  the  interface.  Beyond  this  the  theory  has  serious 
limitations.  Later  developments  of  the  penetration  and  the  boundary 
layer  theories  have  provided  more  realistic  models. 


For  convenience  and  by  usage  it  has  been  customary  to 
express  the  mass  transfer  rate  by  the  equation 


N  =  K(Cb  -  Cs) 


where,  N  =  Rate, 

K  =  Mass  Transfer  Coefficient, 
C-)-,  =  Bulk  concentration, 

Cg  =  Interface  concentration, 


(ii.  i) 

-1  -P 

gmmoles • hr  •  cm 

cm- hr"! 

gmmoles •ml'1  of  solute 
gmmoles *ml"1  of  solute 


This  equation  itself  is  a  mathematical  definition  of  the  mass 
transfer  coefficient  and  this  coefficient  remains  to  be  related, 
theoretically  and  experimentally,  to  the  basic  properties  of  the 
fluid  and  the  flow  conditions. 


Usually  the  experimental  data  on  mass  transfer  coefficient 
are  related  with  the  apparent  variables  of  a  system  by  a  semi- 


' 
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empirical  equation  of  the  form 

ki  /asN.i?f  if i 

Dl  "  t  I  (SI>1 

where,  Dq  =  Diffusivity  in  the  liquid  phase 
1  =  A  characteristic  length 

dg  =  Agitator  diameter 


(II.  2) 


The  equation  can  he  derived  by  dimensional  analysis. 
The  constants  a,  b  and  c  of  the  equation  are  determined  from 
experimental  results.  Though  applicable  only  under  limited 
conditions,  the  equation  has  been  widely  used  to  predict  the 
mass  transfer  coefficient. 

II-  C.  INTRA-PARTICLE  DIFFUSION  AND  SURFACE  REACTION 


From  the  particle  surface  reactants  must  diffuse  into 
the  catalyst  particle  through  the  pore  structure.  The  pores  are 
substantially  interconnected  channels  of  irregular  shapes  and 
varying  cross-sections.  Any  mathematical  treatment  of  intra¬ 
particle  diffusion  remains  largely  emperical  because  the  pore 
geometry  of  most  catalysts  cannot  be  adequately  defined.  A 
simple  approach  has  been  to  represent  the  complex  diffusion  process 
by  a  single  effective  diffusion  coefficient,  D,  which  satisfies  the 


equation 
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(II-  3) 


at  any  location  in  the  particle. 


where,  N 


C 


x 

D 


Flux  of  the  flow  of  diffusing  species, 

gmmoles*min“]-  cm 

Local  concentration  of  the  diffusing 
species  in  the  porous  solid,  gmmoles •  cm”^ 

Distance  parameter,  cms 

Effective  diffusivity  in  the  solid  phase, 

2  .  -1 
cm*  mm 


It  should  be  noted  that  the  effective  diffusivity  is 
based  on  the  unit  geometric  area  of  the  solid  particle. 

As  the  reactants  reach  the  interface,  chemical  reaction 
takes  place.  Many  theories  have  been  developed  in  an  attempt  to 
explain  the  nature  of  surface  reaction.  Very  detailed  mechanisms 
have  been  proposed.  The  most  commonly  accepted  theory  postulates 
that  the  reaction  takes  place  on  active  sites  on  the  surface  of 
the  catalyst.  The  reaction  is  therefore  accompanied  by  adsorption 
and  desorption  processes.  The  rate  expressions  resulting  from 
these  mechanisms  involve  many  arbitary  constants.  Many  a  times  a 
number  of  alternative  mechanisms  fit  the  data  equally  well  and  the 
choice  of  the  correct  reaction  scheme  is  not  obvious.  For  this 
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reason  some  authors  have  advocated  the  use  of  simple  equations 
of  the  form 


Rate  =  kCn 
or 

aC 

Rate  =  - 

1+bC 

to  correlate  rate  data.  This  type  of  rate  expression  simplifies 
the  mathematical  treatment. 

The  simultaneous  mass  transfer  and  chemical  reaction 
establishes  a  concentration  gradient  within  the  catalyst  particle. 
The  interior  surfaces  are  then  exposed  to  reactant  concentrations 
lower  than  the  ambient  fluid  and  the  overall  reaction  rate  through¬ 
out  the  catalyst  particle  under  isothermal  conditions  is  less  than 
it  would  be  if  there  were  no  intra-particle  diffusion  limitations. 

In  this  sense  the  internal  surface  is  not  effectively  utilized. 

An  'Effectiveness  Factor',^,  can  be  defined  as  the  ratio  of  the 
actual  reaction  rate  to  that  which  would  occur  if  all  the  internal 
surface  were  exposed  to  reactant  of  the  same  concentration  and 
temperature  as  that  existing  at  the  outside  surface  of  the  particle. 

Rate  observed 


Rate  in  absence  of  internal  diffusion 


(II. 4) 


10 


The  effectiveness  factor  describes  the  rate  of  reaction  in  terms 
of  the  values  of  the  concentration  and  temperature  at  the  external 
surface. 


For  simple  cases  the  effectiveness  factor  can  be  quantita¬ 
tively  expressed  as  a  function  of  the  catalyst  size  and  intrinsic 
reaction  parameters;  e.g.,  for  spherical  particles  with  no  internal 
temperature  gradients  and  first  order  reaction  (See  Appendix  A) 


1 


— 2 —  (wR-coth(wR)  -  l) 
w2R2 


where, 


k 

D 


(II. 5) 


This  equation  has  important  practical  implications 

:  If  wR  is  small  -  that  is  small  k  and/or  R,  and/or  .large 
D  -  approaches  unity.  All  the  active  surface  is 
exposed  to  the  same  reactant  concentration;  internal 
diffusional  effects  are  insignificant. 

:  If  wR  is  large  -  in  a  system  with  .low  diffusivity  and/or 

fast  reaction  and/or  large  particle  diameter  -  the 
effectiveness,"^  ,  becomes  small,  i.e.,  the  internal 
area  is  not  effectively  utilized. 


1.1 


Thus  very  active  catalysts  (large  k)  tend  to  have  low 
effectiveness  factors.  In  this  case,  maximum  use  is  made  of  a 
catalyst  by  employing  the  smallest  practical  particle  size. 


CHAPTER  III 


LITERATURE  REVIEW 

Ion  exchange  resins  have  been  found  to  catalyse  numerous 
reactions  (40).  Ester  hydrolysis  using  a  cation  exchange  resin 
has  been  subject  to  considerable  investigation.  The  fact  that 
an  acidic  ion  exchange  resin  can  catalyse  the  hydrolysis  of  an 
ester  was  observed  by  Sussman  (55)  and.  Thomas  and  Davies  (57). 

Hammett  et  al  (4,29)  studied  the  hydrolysis  of  ten  esters, 
including  Methyl  and  Ethyl  Acetates,  in  70$  aqueous  acetone  at 
25°  C.  They  used  IR-.120  as  the  catalyst.  Their  batch  experiments 
in  test  tubes  indicated  that 

The  resin  catalyst  was  less  effective  than  using 
a  mineral  acid  as  catalyst,  e.g.  for  ethyl  acetate 

kdesl?)  =  0  226 

k(HCl)  J 

*  Diffusion  rates  were  fast: 

(a)  reaction  rate  was  insensitive  to  the  degree 
of  agitation, 

(b)  particle  diameter  of  the  resin  did  not 
seem  to  have  an  effect  on  the  reaction  rate 

i.e.,  the  chemical  reaction  rate  was  the  limiting  step. 


3V£li 
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*  The  reaction  rate  was  proportional  to  the  resin  mass. 

*  The  rate  constant  k  decreases  with  increase  in 
molecular  weight  of  the  ester.  For  ethyl  acetate 

k  =  1,5  x  10-5  litre* gmmole"'^  sec"-1-  at  25°  C 

In  variance  with  these  results,  Davies  and  Thomas  (l6) 
observed  that  a  cation  exchange  resin  was  more  efficient  than  a 
homogeneous  system  for  the  hydrolysis  of  various  esters  including 
ethyl  acetate.  They  attributed  the  discrepancy  with  Hammett  (4) 
to  the  presence  of  acetone  in  the  latter's  reaction  mixture.  They 
reported  the  following  rate  constants  for  ethyl  acetate,  using 
IR-100  as  the  catalyst. 


Temperature  °C  15 

25 

35 

45 

55 

k  x  icfi  gm“^*sec”^-  3*5 

9.4 

22.7 

49.2 

97.8 

and  E  =  15,540  cal*gmmole 


k(resin) 

k(HCl) 


2.3  at  25°  C 


They  confirmed  that  the  particle  diameter  had  no  effect  on  the 
reaction  rates. 


Bernhard  and  Hammett  (5)  studied  the  hydrolysis  of  methyl 
and  ethyl  acetates  in  aqueous  solutions  using  ion  exchange  resins. 
They  showed  that  reaction  rates  depend  on  the  resin  type.  While 


,3  >L  1  fr  *s 
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a  highly  crosslinked  resin  like  Dowex  -  50  x  8  was  less  effective 
than  HC1  as  the  catalyst,  the  .lightly  crosslinked  resin,  DVB  -  4, 
was  more  effective.  They  reported  the  following  rate  constants 
for  ethyl  acetate. 

k  x  105  litre •gmmole-1  sec"-1-  E  cal*gnmiole"^ 


Temperature  °C 

10 

20 

30 

Catalyst 

HC1  (0.466N) 

2.32 

6.60 

17.20 

17.20 

DVB  -  4 

3.20 

9.21 

24.00 

17.40 

DOWEX  -  50 

2.09 

5.99 

16.00 

17.50 

Hammett  et  al  (6,50)  made  a  detailed  investigation  of 
the  specificity  of  the  reactions  and  showed  that  the  catalytic 
specificity  of  a  resin  catalyst  can  he  altered  by  fundamentally 
two  different  ways:  varying  the  degree  of  crosslinking  and 
changing  the  nature  of  the  functional  group. 

Levesque  and  Craig  (42)  while  studying  the  esterification 
of  some  alcohols  claimed  that  the  reaction  rates  were  proportional 
to  the  resin  surface  and  not  the  mass,  thus  indicating  the  presence 
of  strong  internal  diffusion. 

Saletan  and  White  (52)  studied  a  continuous  process  for 
the  esterification  of  ethanol  with  acetic  acid  by  passing  the 
reactants  through  beds  of  a  cation  exchanger,  Dowex  -  50,  of 
varying  particle  size  at  different  flow  rates  and  temperatures. 
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Their  results  showed  that  the  reaction  rate  was  nearly  independent 
of  the  liquid  flow  rate  even  in  the  .laminar  region.  Hence  liquid 
film  diffusion  effects  were  insignificant.  Their  second-order  rate 
expression  gave  the  average  volumetric  rate  of  reaction  in  a 
spherical  catalyst  particle  as  the  product  of  the  reaction  rate 


at  the  resin  surface  and  the  effectiveness  factor. 
Toy  equation  (II.5): 


,  as  given 


(wR-coth(wR)  -  l)  (II. 5) 

They  suggested  that  the  above  expression  applied  to  close 
approximation  to  second-order  reactions  in  spherical  particles. 

In  their  experiments  the  particle  diameter  was  varied  from  0.009 
to  0.05  cms  and  the  temperature  from  30  to  70  degrees  C;  the 
calculated  effectiveness  factors  varied  from  0.4  to  1.0.  The  results 
show  that  unlike  other  porous  catalysts,  the  diffusivity,  D,  was 
an  exponential  function  of  temperature.  The  rate  of  increase  of 
the  specific  reaction  velocity  constant  with  temperature  was, 
however,  greater  than  the  rate  of  increase  of  diffusivity  with 
temperature.  Therefore  at  high  temperatures  the  internal  diffusion 
was  the  rate  limiting  step,  particularly  for  the  larger  particle 
size  of  the  catalyst. 

In  a  similar  study  Barker  and  White  (2)  measured  conversions 
for  the  alcoholysis  reaction  of  ethanol  and  n-butyl- acetate.  They 
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confirmed,  in  general,  the  observations  of  the  previous  investigators. 
However,  three  entirely  different  mechanisms  correlated  their  data 
within  the  experimental  error. 

Smith  and  Amundson  (54)  considered  the  problem  of  intra¬ 
particle  diffusion  for  catalytic,  simple  reversible  reaction,  first- 
order  in  each  direction.  In  the  theoretical  study  they  derived 
equations  to  relate  the  conversion  with  the  apparent  variables  of 
the  system.  Taking  into  account  the  effects  of  mass  transfer  and 
chemical  reaction  rate  in  the  resin  phase,  they  obtained  expressions 
for  conversion  in  CFSTR,  fixed  bed  and  batch  reactors.  The  mathematical 
models  were  experimentally  verified  for  the  CFSTR  and  batch  reactor. 

The  reaction  employed  was  the  hydrolysis  of  ethyl  formate  catalysed 
by  the  Dowex  -  50  resins.  The  reaction  was  carried  out  at  25  degrees 
C  and  600  r.p.m.  in  a  one-liter  glass  reactor.  The  intense  agitation 
eliminated  the  possible  effects  of  external  film  diffusion  and  non¬ 
ideality  of  mixing. 

For  the  irreversible  reaction  of  pseudo-first  order,  they 
verified  the  linear  dependence  of  C^/Aq,  inlet  ester  concentration 
divided  by  the  acid  concentration  in  outlet,  upon  the  flow  rate, 
q  (See  equation  (lV.l8)  ),  Four  particle  sizes  from  0.005  to  0.08 
cms  were  used.  Their  data  indicate  that  the  effectiveness  factor, Yj  , 
was  nearly  unity  for  the  smaller  particles  and  it  decreased  with 
an  increase  in  the  particle  size.  For  the  hydrolysis  of  ethyl 


formate  they  reported 


■ 


•  q-  -‘  '•  ■  - 


?  ■ 


■ 


k  =  0.962 

D  =  5.0  x  10“ 5 


min 

cm“*min~^  at  25°  C 


Engel  and  Hougen  (l8)  carried  out  hydrolysis  of  amyl 
acetate  in  a  batch  reactor  using  three  different  sizes  of  acidic 
ion  exchange  resins.  They  resolved  the  total  resistance  to  the 
overall  reaction  into  the  liquid  phase  resistance  and  intra¬ 
particle  resistance.  Their  data  were  correlated  by  an  equation 
of  the  form 


Overall  Resistance  = 


A'  (%e)a(HSo)b 


wher  e , 

A’,  B’,  C’,c,  a  and  b  are  constants 

d  =  Particle  size 
P 


and 

Reaction  rate 


Ester  concentration  in_the  bulk  phase 

Overall  Resistance 


They  observed  that  only  for  the  particle  size  of  0.02 6 
mms  was  the  liquid  film  diffusion  significant  and  for  larger 
particles  the  constant  a‘  was  negligible.  They  defined  a 
contactor  efficiency,  0,  as  a  measure  of  the  effectiveness  of 
agitation.  As  f>  approaches  unity,  the  external  mass  transfer 
resistance  decreases  and  the  internal  diffusion  and  chemical 


reaction  became  the  rate  controlling  step. 
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For  a  heterogeneous  reaction  in  a  C’FSTR  the  effect  of 
agitation  upon  reaction  rate  could,  he  complicated,  by  the  occurrence 
of  nonidea.1  flow,  i.e.,  deviation  from  the  assumption  of  complete 
mixing.  This  aspect  has  not  been  studied,  previously. 

The  first  systematic  attempt  to  develop  a  theory  of 
residence  time  distributions  was  made  by  Gilliland  and  Mason  (22,  23) 
and  by  Dankwerts  (.13),  who  introduced  some  quantitative  measures 
relating  to  various  features  of  flow  systems.  Since  then  consider¬ 
able  work  has  been  done  on  this  aspect  (ll,  l4,  15,  46,  60) .  Leven- 
spie.1  (4l)  has  summarized  the  theories  of  residence  time  distribu¬ 
tions,  the  methods  of  interpretation  of  the  actual  experimental 
flow  patterns  and  the  approach  of  postulating  flow  models  for 
prediction  of  the  reactor  performance.  Cholette  and  Cloutier  (.10) 
considered  various  types  of  models  for  a  CFSTR,  derived  related 
theoretical  equations  and  experimentally  verified  one  of  the 
simpler  models. 

Not  many  experimental  observations  have  been  reported  on 
the  measurement  of  actual  residence  time  distributions  for  CFSTR 
measured  as  the  response  to  a  known  stimulus.  However,  the 
techniques  of  measuring  the  residence  time  distributions  and  the 
relative  advantages  of  various  methods  are  extensively  treated  by 
Levenspiel  (4.l)  and  Naor  and  Shinnar  (46).  Cholette  and  Cloutier 
(.10  )  explored  the  nonideal  behaviour  of  a  30  inches  dia.  x  30 
inches  height  stirred  tank  and  found  that  their  data  was  best 


. 
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described  by  combination  of  a  backmix  and  a  dead  water  region  with 
a  portion  of  the  feed  short-circuiting  the  vessel.  They  measured 
the  response  of  the  stirred  tank  to  a  step  demand  of  water  in  a 
0.05  N  NaCl  solution.  The  NaC.1  concentrations  were  measured  at 
regular  time  intervals  by  volumetric  titration.  For  their  experi¬ 
mental  conditions  nearly  ideal  flow  prevailed  at  agitator  speeds 
of  more  than  200  rpm.  No  other  reference  can  be  cited  in  support 
of  the  above  work.  In  contrast.  Wolf  and  Manning  (59)  measured 
the  residence  time  distributions  of  a  stirred  tank  over  the  range 
of  0-500  rpm  and  the  data  indicate  that  it  changes  very  little  with 
the  stirrer  speed.  Also,  Worrell  and  Engles ton  (6.l)  noticed 
insignificant  difference  between  the  residence  time  distributions 
for  no  mixing  and  perfect  mixing;  they  measured  the  response  of  a 
pulse  injection  of  10  °]o  KC1  solution  in  the  flow  of  distilled  water. 
Other  investigators  (7,38)  have  shown  that  homogeneity  can  be 
achieved  in  as  short  a  time  as  one  second  at  20-50  rpm  in  a  .laboratory 
sized  stirred  tank. 

In  the  field  of  external  film  diffusion,  Harriott  (27,  28), 
Hyman  (36)  and  Miller  (44)  have  given  critical  reviews  of  the 
theoretical  studies  on  solid-liquid  mass  transfer  in  agitated 
vessels.  Miller  (44)  and  Marangozis  and  Johnson  (43)  have  collected 
available  experimental  data  on  the  solid-liquid  systems  and  compared 
various  correlations  proposed  in  the  literature  (2,  9,  21,  25,  28, 

35,  37,  45,  62). 
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Most  of  the  measurements  have  been  made  on  the  processes 
involving  the  dissolution  of  a  solute  which  was  in  the  form  of 
freely  suspended  particles.  Harriott  (27)  and  Calderbank  and 
Moo-Young  (9)  measured  the  mass  transfer  coefficients  with  ion 
exchangers.  The  effect  of  vessel  geometry  and  the  properties  of 
the  solid-liquid  systems  used  have  been  studied  and  the  experiments 
cover  a  wide  range  of  variables.  The  data  have  been  mostly  correlated 
by  the  Gilliland  -  Sherwood  equation 


C(N  )a 
K  Re' 


(N_  ) 

'  Sc' 


b 


(II. 2) 


In  most  of  the  correlations  the  exponent  b  has  been 
assumed  to  be  equal  to  0.5  to  concur  with  the  penetration  theory. 
However,  Marangozis  and  Johnson  (43)  have  shown  that  most  of  the 
data  can  also  be  fitted  by  b  =  l/3,  as  predicted  by  the 
boundary  layer  theory.  Both  the  theories  seem  to  apply. 

The  exponent  a  varies  from  0.4  to  0.9  among  different 
investigators.  Harriott  (27)  has  proposed  an  alternate  approach. 

He  argues  that  the  use  of  stirrer  diameter  as  the  characteristic 
length  in  the  Reynold  and  Sherwood  numbers  is  basically  erroneous 
and  consequently  these  correlations  have  no  fundamental  significance. 
He  justifies  the  use  of  particle  diameter  instead.  Mass  transfer 
to  and  from  suspended  solids  depends  primarily  on  slip  velocity, 
with  the  terminal  velocity  yielding  the  maximum  value  of  the  mass 


- 
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transfer  coefficient,  which  can  be  calculated  by  using  Fr os sling’s 
semi-theoretical  equation 

TT  =  2  +  °-55(NEe)l/2(NSc)l/3 

The  actual  mass  transfer  coefficient  for  a  given  system  can  then 
be  calculated  from  his  graphs  for  K/K*  versus  power  input. 

Presently  available  correlations  are  uncertain  because  of 
the  discrepancies  in  the  values  of  the  constants,  the  accompanying 
scatter  of  the  data  and  the  wide  differences  in  the  reported  effects 
of  the  variables  like  particle  size  (3,  l8,  28,  30-35?  58),  stirrer 
speed  (3?  24,  35)?  stirrer  location  (28,  37)  and  baffles  etc  (3?  19? 

21,  30,  36,  37). 

Following  the  classical  work  of  Thiele  (56)  numerous  articles 
have  appeared  on  the  subject  of  internal  diffusion  and  chemical 
reaction  in  porous  catalysts.  From  this  vast  literature  has  evolved 
a  mathematical  treatment  which  has  now  become  a  standard  approach 
to  the  problem.  Excellent  books  have  also  appeared  on  the  subject 
(47?  53)*  As  a  result,  no  detailed  review  of  the  literature  was 
undertaken.  While  on  the  theoretical  side  massive  work  is  being 
done  to  explain  the  phenomenon  of  catalysis  (.1,  39)  and  in  deriving 
rate  expressions  for  complex  reaction  mechanisms  and  pore  geometries 
(.12,  20,  39?  47)  the  most  common  approach  on  the  engineering  side 
has  been  to  use  an  average  effective  diffusivity  and  a  simple  rate 
expression. 
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CHAPTER  IV 


REACTION  SYSTEM  MODEL 

When  a  solid-liquid  catalytic  reaction  is  carried  out  in 
a  CFSTR  the  intensity  of  turbulence  in  the  flow  field  affects 
simultaneously  the  processes  of  transport  and  mixing,  the  former 
being  the  transfer  of  matter  under  a  concentration  gradient  and 
the  latter  the  rate  of  progress  of  the  system  toward  the  equili¬ 
brium  state  of  complete  randomness  at  the  molecular  level.  Implicit 
with  the  concept  of  an  ideal  CFSTR  is  the  assumption  that  the 
contents  of  the  reactor  are  so  thoroughly  mixed  that  they  are  at 
uniform  concentration  and  temperature  throughout.  However,  insuffi¬ 
cient  mixing  may  occur  at  low  impeller  speeds  which  causes  deviations 
from  the  ideal  behaviour  of  the  reactor.  The  concentrations  are 
then  not  uniform  in  the  reaction  mixture  and  the  observed  rate  of 
the  reaction  differs  from  that  obtainable  in  an  ideal  CFSTR. 

Many  models  have  been  proposed  to  characterize  nonideal 
flow  patterns  so  that  quantitative  prediction  of  the  performance 
of  real  equipment  can  be  made.  The  usual  procedure  is  to  obtain 
the  actual  residence  time  distribution  for  the  vessel  and  to  use 
this  information  for  applying  correction  to  the  ideal  flow.  The 
residence  time  distribution  gives  a  measure  of  the  time  individual 
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molecules  stay  in  the  vessel.  A  residence  time  distribution 
function  F(t)  of  a  continuous  flow  system  is  defined  as  the  volume 
fraction  of  the  fluid  at  the  outlet  which  has  resided  in  the  system 
for  a  time  less  than  t.  This  function  F(t)  represents  the  proba¬ 
bility  that  an  element  of  volume  which  has  entered  the  system  at 
t=0 ,  has  left  it  within  a  period  of  time  t;  the  probability  that 
it  will  leave  a  moment  .later  than  t  is  1  -  F(t).  Therefore,  dF(t) 
is  the  volume  fraction  of  the  outgoing  stream  which  has  a  residence 
time  between  t  and  t  +  dt. 

Clearly, 


t  =  0; 


and  average  residence  time 


F(t)=0 

F(t)=l 


X) 


Levenspia.1  (4l)  discusses  other  residence  time  distribution 
functions.  For  idealized  systems  these  functions  can  be  derived 
mathematically.  In  principle,  all  the  information  is  embodied  in 
any  one  of  the  various  functions  and  in  systems  where  steady  state 
conditions  prevail  the  relations  among  them  are  uniquely  determined 
by  mathematical  equations. 


For  an  ideal  backmixed  vessel  it  can  be  shown  that 


-© 


if  0  =  t/t 

F(e)  =  1  -  e 


(IV.  1) 
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Experimental  measurements  of  the  residence  time  distribution  are 
done  by  the  stimulus-response  technique,  that  is  measuring  the 
response  of  the  vessel  to  a  known  input  signal  of  a  tracer.  Most 
commonly  a  step  function  or  a  delta  function  is  used.  It  can  be 
shown  that  for  a  perfectly  backmixed  vessel,  the  response  to  an 
imposed  step  function,  measured  as  the  concentration  time  curve 
at  the  vessel  outlet  is  F(q). 
i.  e. , 

F(Q)  =  =  i  -  e"0  (IV. 2) 

cn(o) 


where, 


CN(©)  =  Tracer  concentration  at  time  © 


This  information  can  be  used  directly  only  for  first  order  homo¬ 
geneous  reactions.  For  nonlinear  reaction  rates  it  is  necessary 
to  hypothesize  what  is  considered  to  be  a  reasonable  model  for 
the  behaviour  of  the  fluid  in  the  vessel  and  then  in  conjunction 
with  the  residence  time  distribution  information  calculate  the 
conversions.  The  flow  model  approach  consists  in  visualizing  the 
vessel  with  imperfect  mixing  as  made  up  of  two  or  more  flow 
regions  such  as 

Plug  flow  regions 
Backmix  regions 
Dead  water  regions 
By-pass  regions,  etc. 
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For  each  combination  there  exists  a  theoretical  response 
curve.  It  can  be  derived  mathematically  in  terms  of  the  para¬ 
meters  which  express  the  magnitude  of  the  constituent  regions  of 
the  reactor  model.  The  theoretical  response  curve  is  then 
matched  as  closely  as  possible  to  the  experimental  response 
curve.  This  yields  the  values  of  the  various  parameters. 

The  problem  of  nonideality  in  a  heterogeneous  system 
could  be  complicated.  In  a  backmix  reactor  with  no  through  passage 
of  solids  nonideal  behaviour  can  exist  due  to  nonuniform  distribu¬ 
tion  of  solids  in  the  reactor,  in  addition  to  the  nonideal  flow  of 
liquid.  Super ic ally  this  can  be  taken  into  account  by  proposing  a 
flow  model.  A  simple  model  is  proposed  for  the  system  under  study. 
A  fraction,  m,  of  the  feed  is  assumed  to  by-pass  the  solids, 
while  in  the  rest  of  the  reactor  volume  the  solids  are  uniformly 
distributed;  that  is  to  say  that  the  reactor  is  composed  of  a  by¬ 
pass  zone  in  which  there  are  no  solids  and  a  perfectly  backmixed 
zone  which  contains  all  the  catalyst.  There  is  no  conversion  in 
the  by-pass  stream,  while  in  the  other  stream  the  conversion  is 
that  obtainable  in  an  ideal  backmix  reactor.  At  the  outlet,  it  is 
hypothesized,  the  two  streams  mix  and  the  resulting  composition 
can  be  calculated  by  material  balance.  A  pictorial  representation 
of  the  model  is  shown  in  FIGURE  (iV-l) . 

With  the  simple  model  of  by-pass  plus  backmix  flow,  it  is 
possible  to  derive  expressions  which  can  be  used  to  determine  the 


' 


. 


o  p 


26 


FIGURE  IV-  I  *  REACTOR  MODEL 
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parameters  of  the  reaction  system.  In  the  following  analysis  the 
reaction 

Resin 

ch3cooc2h5  +  h2o  _ _  ch3cooh  +  c2h5oh 

takes  place  in  a  reactor  containing  W  gms  of  catalyst  and  V 
liters  of  solution.  The  feed  is  introduced  at  a  rate  of  q  mis 
per  minute  and  has  an  ester  concentration  of  Cj_  gmmoles/ml.  Effluent 
solution  is  withdrawn  at  the  same  rate  hut  with  an  acid  concentration 
of  AQ  and  the  ester  concentration  of  C3.  The  feed  solution  does 
not  contain  any  acid,  therefore,  at  steady  state 

ci  =  Co  +  Ao  (IV. 3) 

According  to  the  model  proposed,  only  a  fraction  m  of 
the  feed  enters  the  zone  of  perfect  mixing,  while  the  other  fraction 
(.1  -  m)  is  assumed  to  be  short-circuited  directly  to  the  outlet. 

The  concentration,  at  steady  state,  in  the  zone  of  perfect  mixing  is 
assumed  to  be  C*  .  The  material  issuing  from  this  zone  is  mixed  with 
that  portion  of  the  feed  which  short-circuits  the  system;  the  mixture 
of  the  two  giving  a  liquid  of  composition  CQ  at  the  outlet. 

A  material  balance  at  the  outlet  gives 


(l  -  mJCq  +  mC* 


(iv. 4) 
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The  reaction  is  considered  to  be  elementary,  irreversible 
and  represented  kinetically  by 


Rate  =  kC 


(IV. 5) 


-1 


C  =  Concentration  of  the  ester  gmmoles*ml 

k  =  Chemical  rate  constant,  min'"-'-;  it  is  characterized  by  the 
catalyst. 


For  a  single  sphere,  a  material  balance  at  steady  state 
over  a  spherical  shell  of  thickness  dr  gives 


(4nr2  D 

dr  r  +  dr 


(4nr2  D  <E) 
dr  r 


—  (4nr?dr*kC)p  =  0 

where, 

D  =  Effective  diffusivity  of  ethyl  acetate  in  the  resin  phase, 

cm?min“-'- 


as  dr 


0,  this  equation  reduces  to 


dfC 

dr^ 


+  1 


ac 

dr 


kC 

D 


=  0 


Substituting 


and 


u  =  rC 


k  2 
D  -  w 


(IV. 6) 
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in  the  above  equation 

—  -  w2u  =  0  (IV.  7) 

dr2 


A  General  solution  of  the  equation  (IV. 7)  is 


u  =  Acosh(wr)  +  Bsinh(wr) 


A  B 

C  =  —  cosh(wr)  +  _  sinh(wr) 

r  r 


Boundary  conditions  for  the  equation  are 

(l)  At  the  centre  of  the  catalyst  particle 


=  0 


(IV. 9) 


(2)  At  the  catalyst  surface,  the  ester  diffusing  through 
the  liquid  film  enters  the  particle. 

D  f^£)  =  K(C*  -  C_)  (IV.  10) 

l dry  r=R  S 

where 

K  =  Mass  transfer  coefficient,  cms-min"1 

Cg  =  Concentration  of  ester  on  the 

particle  surface,  gmmoles*ml  1 

R  =  Average  radius  of  the  catalyst 

particles  as  defined  by  equation  (V.l),  cms 
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The  rate  at  which  the  ester  enters  a  single  sphere  is 


4rtR2  d/— \ 

W)  r=R 


and  the  rate  at  which  it  enters  W  gms  of  spheres  is 


3WD 

RS 


r=R 


(IV.  u) 


where 


s  = 


w 


Volume  of  W  gms  of  wet  resins  ml 


Hence,  this  is  the  rate  at  which  the  main  body  of  the  solution  is 
being  depleted  in  ester  by  diffusion  into  the  particles. 


q(Ct  -  C*)  =  22  /dC\ 

R?  ld r1  r=E 


(IV. 12) 


Eliminating  C*  from  equations  (IV. 12)  and  (IV. 10) 


Ci  -  Cs 


3W  D 

R^mq 


+  R 

K 


'dCA 

,drj  r=R 


(iv.  13) 


The  differential  equation  (IV.6),  in  conjunction  with  the 
boundary  conditions  (IV. 9)  and  (IV. 13)  yields  the  solution  in 
the  following  form  (See  Appendix  A). 


Ci 

A 


7-  =  <1 


O 


1  + 

Q, 


R£_ 

3KW 


+  1 
m 


(iv. 14) 
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where 


(iv. 15) 


where 


w 


D 


It  is  possible,  in  principle,  to  evaluate  all  the  four 


parameters  K,  m,  k  and.  D  from  more  than  four  experimental  values 
C  • 

of  at  varying  R,  W  and  q  and  constant  N  and  T,  by  using  equation 

Ao 

(IV.14) .  Mathematically,  however,  it  would  be  a  staggering  task 
due  to  the  complexity  of  the  equation.  A  simpler  and  more  reliable 
method  of  evaluating  the  parameters  consists  of  two  steps:  First 
determine  k  and  D  and  then  evaluate  K  and  m  at  desired  conditions. 

The  method  is  elaborated  below. 


When  the  ratio  of  the  inlet  ester  concentration  to  the 


exit  acid  concentration,  _1,  in  an  experimental  run  is  plotted 

Ao 

versus  the  flow  rate,  q,  equation  (lV.l4)  gives  a  straight  line  with 


Slope 


1  +  Ri 


(iv. l 6) 


Q  3KW 


and 


Intercept  =  — 

m 


(IV. 17) 


As  can  be  seen  from  the  equation  (IV.,15),  for  a  given 


quantity  of  a  catalyst,  Q  is  a  function  of  only  temperature 
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and  the  particle  diameter  and  is  independent  of  the  reactor 
geometry  and  the  flow  conditions.  To  evaluate  the  constants 
k  &  D,  two  conditions  in  the  reactor  operation  need  to  be  satisfied: 

(a)  The  reactor  behaves  as  an  ideal  CFSTR,  and 

(b)  The  mass  transfer  in  the  external  film  is  very  rapid. 

If  one  or  both  of  the  above  conditions  are  not  met,  the 
reaction  rate  would  be  considerably  affected  by  changes  in  the 
intensity  of  turbulence  in  the  fluid  field.  A  simple  method  of 
ensuring  that  the  above  conditions  are  met  is  based  on  the  expectation 
that  at  high  agitation  rates  the  increased  turbulence  must  eventually 
became  adequate  to  supply  reactant  at  the  rate  at  which  it  would  be 
sufficient  for  consumption  at  the  catalyst  surface  and  also  to  give 
complete  mixing  of  the  reactor  contents  so  that  an  ideal  backmix 
flow  is  approximated.  When  the  reaction  rate  is  constant  as  agitation 
is  varied,  the  concentration  in  the  bulk  phase  is  assumed  to  be  uniform 
and  it  is  reported  that  no  nonideality  in  flow  exists,  i.e.,  m  =  1; 
it  is  also  assumed  to  be  the  same  as  at  the  catalyst  surface  and  it 
is  reported  that  no  film  mass  transfer  resistance  exists,  i.e.,  K  =  0°  . 

When  the  above  conditions  are  met,  the  equation  (IV.14) 
simplifies  to 

Ci 
Ad 


1  + 


£ 

Q 


(iv. 18) 
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Ci 

The  plot  of  —  versus  the  flow  rate,  q,  should,  give 
Ao 

a  straight  line  with 

Slope  =  i. 

Q 

Intercept  =  .1 

From  a  given  set  of  the  experimental  values  of  Q  ,  the 
parameters  k  &  D  can  then  be  calculated  from  the  equation 

Q,  =  (wR.coth(wR)  -  l)  (IV.  15) 

R^5 

By  varying  the  temperature  in  the  above  procedure  the 
rate  constants  can  be  calculated  as  a  function  of  temperature. 

The  experiments  can  now  be  performed  under  any  desired 
conditions  and,  as  stated  above,  the  equation  (lV.,l4)  plotted  to 
determine  the  slope  and  the  intercept.  Q  calculated  from 
equation  (IV.,15)  is  substituted  in  the  expression 

Slope  =  I  +  _M  (IV.  l6) 

Q  3KW 

to  calculate  the  mass  transfer  coefficient,  K.  The 

Intercept  =  —  (IV.17) 

m 

directly  gives  the  parameter  m.  These  quantities  can  be  determined 
in  relation  to  any  important  variable  of  the  system,  e.g.,  stirrer 
speed,  particle  diameter  etc. 
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The  values  of  K  &  m  thus  obtained  can  be  compared  with  those 
evaluated  in  a  different  manner.  The  coefficient  of  mass  transfer 
can  be  calculated,  as  stated  previously,  by  a  suitable  equation 
of  the  form 

Ksh  =  c(%e)a(NSc)  (II.  2) 

The  parameter  m  can  be  evaluated  directly  by  the 
stimulus-response  technique.  It  has  been  shown  (4l)  that,  for 
the  proposed  model,  the  theoretical  response  curve  to  a  step 
function  in  the  inlet  concentration  is 

F(9)  =  C-.(.Q)  =  1  -  me‘mG  (IV.19) 

CN(O) 

The  function  F(©)  is  obtained  experimentally. 

Once  the  adequacy  of  the  proposed  model  is  established 
for  an  actual  reactor,  equation  (lV.,l4)  can  be  used  for  various 
purposes,  e.g. 

(1)  To  predict  the  actual  conversion  in  the  reactor 
by  using  known  K,  m,  k  and  D  for  given  R, 

W  and  q. 

(2)  To  determine  the  rate  parameters  k  and  D  from 
experimental  observations  at  any  stirrer  speed  if 
values  of  K  and  m  under  the  experimental 
conditions  are  available  by  independent  means. 


CHAPTER  V 


EXPERIMENTAL 


The  experimental  program  was  planned  to 


(1)  establish  the  extent  of  nonideality  in  the  reactor  operation 
at  various  stirrer  speeds, 

(2)  investigate  the  significance  of  liquid-film  mass  transfer 
resistance  on  the  overall  rate  process, 

(3)  evaluate  the  rate  constants  k  and  D  and  if  the  above 
effects  are  significant  to  test  the  validity  of  the  equation 
(IV. 14): 


Ci 

Ao 


R?  +  1 

3KW  m 


(IV. 14) 


V-A  EXPERIMENTAL  EQUIPMENT 

The  experimental  set  up  was  designed  so  that  known  quantities 
of  ethyl  acetate  and  water  could  be  introduced  into  a  CFSTR  in  which 
a  solid  catalyst  was  dispersed  and  the  product  stream  could  be  with¬ 
drawn  from  the  reactor  at  the  same  rate,  retaining  the  catalyst 
inside  the  reactor  by  employing  a  filter. 

The  equipment  consisted  principally  of  feed  reservoirs,  feed¬ 
pumps,  rotameters,  associated  tubings,  a  glass  reactor  placed  in  a 
constant-temperature  water  bath  and  a  proportioning  pump.  A  schematic 
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FIGURE  V- 1  *  SCHEMATIC  DIAGRAM  OF  THE  EXPERIMENTAL  SET-UP 
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diagram  of  the  experimental  set  up  is  shown  in  FIGURE  (Vr.l) 

Introduction_of  feedj_  A  small  centrifugal  pump  supplied  distilled 
water  frcm  a  12-gallon  glass  bottle  to  a  small  glass  bottle  placed 
about  four  feet  above  the  reactor.  The  overhead  bottle  served  as  a 
constant  head  tank  from  which  the  overflow  returned  to  the  feed 
bottle,  and  a  regulated  stream  of  water  flowed  through  a  calibrated 
rotameter  into  the  reactor.  The  tubing  employed  was  inch  diameter, 
stainless  steel. 

A  similar  set  up  was  provided  for  the  ester  stream.  However, 
the  reservoirs  were  smaller  in  size.  A  small  positive  displacement 
pump  served  as  the  feed  pump.  The  ester  dissolves  most  of  the 
commonly  available  rubbers;  therefore,  connections  using  rubber 
stoppers  and  plastic  tubes  were  avoided.  Special  glass-to- steel 
connections  were  used.  Stainless  steel  tubing  of  l/8  inch  diameter 
was  used  for  this  stream. 


Reactor:  The  reactor  vessel  was  constructed  from  a  6  x  12  inch 

cylindrical  pyrex  glass  jar.  The  jar  was  modified  to  give  a  vessel 
with  an  overall  length  of  9  2  inches  and  was  provided  with  a  3/^  inch 
wide  flange.  The  jar  could  then  be  held  between  two  metal  flanges. 
Rubber  gaskets  were  employed  to  provide  a  seal.  The  metal  cover  of 
the  reactor  served  as  the  compliment  ary  part  to  a  stainless  steel 
ring,  6-J-  inches  i.d.  and  1~  inches  wide,  which  served  as  the  lower 
flange.  The  top  cover  was  made  of  a  J  inch  thick  stainless  steel 


' 
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plate,  9  inches  in  diameter  and  it  accomodated  suitable  connections 
for  inlet  and  outlet  stream  lines  and  a  centrally  placed  stirrer 
shaft.  Four  symmetrically  located  stainless  steel  baffles  9ir  inches 
long  and  inch  wide  were  attached  to  the  top  cover  such  that  when 
it  covered  the  glass  reactor  the  clearance  between  the  glass  wall 
and  the  baffles  was  about  l/l6  inch.  In  addition,  the  top  cover 
supported  the  stirrer  assembly  at  a  fixed  position  so  that  the 
stirrer  shaft  was  centrally  positioned  and  accurately  aligned  with 
the  bushing  and  that  the  occasional  removal  of  the  top  cover  did 
not  disturb  the  alignment. 


The  geometrical  features  of  the  reactor  were  chosen  as 
being  typical  of  previous  practice  (36).  The  dimensional  and  other 
details  are  given  in  FIGURE  (V-2)  and  Appendix  B.  Some  important 
dimensions  were: 


I.d.  of  the  reactor 
Liquid  depth 

Volume  of  reaction  mixture 
Impeller  diameter 

Impeller  position 


inches 

^  inches -equal  to  the 
reactor  i.d. 

2  litres. 

13/^  inches  -  equal  to 
l/3  to  4.  of  the  reactor  i.d. 
1  3/^  inches  from  the 
reactor  bottom  -  at  about 
l/3  the  liquid  depth 
■§■  inch  -  equal  to  l/lO  to 
l/l2  the  reactor  i.d. 


Baffle  width 


. 


■*.  noJti  ■ 


LEGEND 

©  BAFFLES  ,  NO  4 
©  ESTER  FEED  LINE,  1/8" 
©  THERMOWELL 
©  WATER  FEED  LINE,  1/4" 
©EFFLUENT  LINE,  1/4" 
©  FLANGE  RING 
®  FILTER 

®  BOLTS  ,  3/8  ,  NO  6 
©  GLASS  REACTOR 


FIGURE  V-2  *  REACTOR  DETAILS 
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Two  types  of  stirrers  were  employed.  An  air  driven  stirrer 
was  used  for  high  impeller  speeds.  It  was  supported  on  a  frame  of 
aluminum  T- section.  The  frame  was  bolted  to  the  reactor  top.  This 
stirrer  was,  however,  unsuitable  for  experiments  at  .low  r.p.m.  for 
it  provided  very  little  torque  at  these  speeds.  For  this  purpose  a 
stirrer  driven  by  an  enclosed  l/50  h.p.  electric  motor,  provided 
with  a  shaft  of  36:1  reduction  and  which  could  deliver  a  torque  of 
9  inch-pound  at  speeds  from  20  to  .175  r.p.m.,  was  used.  The  motor 
was  supported  by  a  cast  aluminum  control  housing  which  in  turn  was 
clamped  to  a  support  rod  mounted  on  the  reactor  top.  The  agitator 
was  a  three  blade,  1  3/4  inches  diameter,  316  stainless  steel  turbine 
impeller. 


A  thermometer  well  made  of  5/16  inch  stainless  steel  tube 
sealed  at  the  lower  end  was  welded  to  the  reacter  top  to  project 
about  2  inches  below  the  liquid  .level  in  the  reactor.  A  mercury 
thermometer  was  placed  inside  the  tube  to  measure  the  reaction  temp¬ 
erature.  The  feed  tubes  terminated  approximately  2  inches  above  the 
liquid  level. 

Product  removal:  A  small  sintered  brass  filter  with  pores  of 

.12-50  micron  size  was  connected  to  a  ^  inch  stainless  steel  tube 
and  placed  in  the  bulk  of  the  reaction  mixture  approximately  2 
inches  above  the  bottom  and  j?  inch  away  from  the  wall  of  the  reactor. 
At  this  location  the  strong  liquid  currents  from  the  impeller  reduced 
the  possibilities  of  filter -choking.  A  -J-  inch  tube-to  ij- inch  tube 
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stainless  steel  bulk-head,  connector  positioned  on  the  reactor  top 
provided  a  leak-proof  continuation  of  the  exit  stream  from  the  inside 
of  the  reactor  to  a  Milton-Roy  controlled  volume  pump  with  a  delivery 
rate  variable  from  0-400  mls/min.  The  pump  was  placed  3  feet  below 
the  reactor  assembly.  As  a  result,  the  liquid  from  the  reactor  was 
siphoned  under  gravity  to  the  pump  inlet.  The  quantity  of  the 
liquid  pumped  was  controlled  by  manipulating  the  stroke  length  of  the 
plungers.  A  back  pressure  valve  adjusted  for  60  psig  was  installed 
in  the  line  near  the  pump  outlet  to  insure  consistent  performance  of 
the  pump.  These  arrangements  permitted  withdrawal  of  a  known  quantity 
of  liquid,  while  retaining  the  catalyst  inside. 

Temperature  control:  A  1 6  x  12  inch  pyrex  glass  jar  filled  with 

water  served  as  the  constant  temperature  bath.  A  1^00  watt  heater, 
an  electric  stirrer  and  a  bimetallic  sensing  element  in  conjunction 
with  an  electronic  controller  were  used  to  maintain  a  constant  water 
temperature  in  the  bath.  A  mercury  thermometer  placed  near  the 
thermoregulator  indicated  the  bath  temperature.  The  set  up  was 
tested  at  125°  F  for  24  hours.  The  temperature  remained  constant 
within  +  0.1  and  the  set  up  was  considered  satisfactory  over  an 
80-120°  F  range. 

The  water  feed  stream  was  preheated  by  passing  a  two  foot 
length  of  the  inlet  line  through  the  water  bath.  This  arrangement 
ensured  that  the  water  feed  was  at  about  the  same  temperature  as 


the  contents  of  the  reactor.  Similar  arrangement  was  not  found 
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necessary  for  the  ester  feed  because  of  its  small  size.  During 
test  runs  excellent  temperature  control  was  attained  and  the 
reactor  contents  were  at  essentially  the  same  temperature  as  the 
water  bath. 

Reactor  support:  A  triangular  frame  of  angle-iron  mounted  on 

the  top  of  the  temperature  bath  jar  held  it  in  a  fixed  position 
and  from  it  the  reactor  was  suspended  in  the  jar.  In  this  position, 
about  8  inches  of  the  reactor  was  immersed  in  the  water  bath. 

The  temperature  bath  and  the  accompanying  assembly  was 
placed  about  3  feet  above  the  ground  level  on  a  bench  to  permit 
convenient  observation  of  the  reactor  operation. 

V-B  REACTANTS  ADD  THE  CATALYST 

V-B-l  Ion  Exchange  Resins: 

Ion  exchange  resins  are  high  molecular  weight  polyacids 
or  polybases  which  are  virtually  insoluble  in  most  aqueous  and 
nonaqueous  media.  In  addition  to  ion  exchange  phenomenon,  their 
utility  encompasses  catalytic  applications.  Since  a  cation 
exchange  resin  owes  its  capacity  for  exchanging  ions  to  an  acidic 
functional  group,  it  may  be  expected  to  promote  reactions  which 
can  be  catalysed  by  conventional  acids. 

Of  the  numerous  types  of  ion  exchange  resins  synthesized, 
the  most  important  types  are  those  prepared  by  suspension  copoly¬ 
merization  of  styrene  plus  variable  amounts  of  divinylbenzene 
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which  is  used  as  a  crosslinking  agent.  The  percentage  of 
divinylbenzene  in  the  backbone  of  the  resin  catalyst  is  commonly 
referred  to  as  the  degree  of  crosslinkage.  These  resins  are 
mechanically  strong  and  chemically  inert.  The  styrene -divinyl¬ 
benzene  network  gives  a  maximum  resistance  to  oxidation,  reduction, 
mechanical  wear  and  breakage  and  is  insoluble  in  common  solvents. 

In  addition,  it  can  be  easily  formed  in  the  shape  of  spheres  or 
beads  of  good  hydraulic  properties.  These  resins  appear  to  offer 
several  advantages: 

•  Simplicity  with  which  they  can  be  removed  from  the  reaction 
medium  by  f iteration. 

•  Elimination  of  the  neutralization  of  the  reaction  mixture. 

•  Acid  catalysis  corrosion  problems  are  reduced. 

•  Cost  reduction:  the  catalyst  can  be  used  repeatedly  without 
regeneration  and  in  many  instances  with  no  loss  of  the 
catalyst  activity. 

They  are  however  unsuitable  for  application  at  temperatures  higher 
than  120°  C. 


The  strongly  acidic  sulfonic  acid  resins,  prepared  by 
treating  the  styrene- divinylbenzene  copolymer  with  hot  concentrated 
sulfuric  acid  are  generally  employed  as  acid  catalysts.  Their 
structure  is 


CH  -  CH, 


SO  H 
3 
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The  important  resins  among  those  commercially  available 
resins  are  Dowex  -  50  and.  Dowex  -  50W,  manufactured  by  the  Dow 
Chemical  Co. ,  and  Anberlyst  -  15  and  IR  -  120,  manufactured  by 
the  Rohm  and  Haas  Co.  Their  properties  and  methods  of  evaluation 
are  available  in  the  corresponding  company  brochures  and  elsewhere. 

(17,  20,  51) 


When  immersed  in  water,  dry  resin  particles  swell  to  an 
extent  determined  by  the  degree  of  cross! inking.  As  it  swells 
into  a  jelly-like  mass,  its  structure  becomes  more  porous  and 
diffusion  within  the  resins  becomes  faster. 


Dowex  50  W  -  x8  was  chosen  as  the  catalyst  because  it  is 
readily  available  in  the  hydrogen  form  and  can  be  obtained  graded 
as  to  size.  The  resins  were  obtained  from  the  Dow  Chemical  Co.  in 
three  size  ranges.  Some  of  the  resin  properties  are: 


Degree  of  Crosslinking 
Active  group 
Ionic  form 
Physical  form 
Size  range 
Sphericity 
Effective  pH  range 
Total  exchange  capacity 


8  i 

Nuclear  sulfonic  H+  SO" 

H+ 

Spheres 

20-50,  50-100,  100-200  mesh 
90  1o 
l-i4 

5  meq/gm  dry  resin 

no.  of  ionic  sites  per  gm  of 


dry  resin. 
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Measurement  of^jurticle  diameter^  Considerable  variation  in 
the  size  of  the  particles  in  the  samples  received  was  observed. 
However  further  attempts  to  obtain  a  sharper  size  separation  met 
little  success  because  very  small  resin  particles  had  a  tendency 
to  adhere  to  relatively  bigger  particles.  The  average  diameter 
for  each  sample  was  therefore  measured  as  the  mean  of  about  one 
hundred  particles. 


by  using  the  equation 


d 


avg 


Xa2 


(v.i) 


This  equation  gives  the  mean  specific  area  diameter,  which  is  the 
diameter  of  a  sphere  having  a  specific  surface  area  equal  to  the 
average  specific  area  of  the  particles.  The  diameter  of  the  wet 
resin  particles  were  measured  by  using  a  microscope. 

About  two  grams  of  wet  resin  was  put  on  a  filter  paper 
to  remove  the  bulk  of  water  from  the  particles.  The  quantity  of 
resin  was  then  symmetrically  divided  on  the  filter  paper  into  eight 
equal  parts.  One  of  the  parts  was  taken  as  a  representative 
sample.  The  particles  were  spread  on  a  glass  slide  and  focused 
under  a  microscope.  The  micrometer  consisted  of  a  movable  table 
with  micrometer  motion  controls  and  an  eye-piece  with  a  slide  wire 
and  a  scale.  To  avoid  any  biased  selection,  only  those  particles 
which  crossed  or  touched  the  horizontal  reference  line  were  taken 
for  measurement.  The  particles  were  essentially  spherical.  The 


' 
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diameter  measured  was  the  horizontal  distance  travelled  by  the  slide 
wire  between  two  diametrically  opposite  vertical  tangents  to  the 
projected  circle.  For  each  sample,  about  one  hundred  particle 
diameters  were  measured.  The  results  are  given  below 


mesh  size 

20  .,50 
50,100 
.100,200 


Average  diameter,  cms. 

0.06090 

0.02536 

O.OO566 


The  details  of  the  measurements  are  given  in  TABLE  (B-l)  of 
Appendix  B. 


Measurement  of_r_esi.n_densityj_  About  ten  gms  of  an  oven- dry 
sample  of  resin  was  accurately  weighed  and  put  in  a  long  graduated 
cylinder.  Water  was  then  slowly  added  for  the  rehydration  of  the 
resin.  After  twenty- four  hours,  the  excess  water  was  drained  off 
the  resins  were  transferred  to  a  filter  paper  and  then  quickly  into 
a  graduated  cylinder  in  which  a  known  volume  of  water  was  present. 

The  resins  were  allowed  to  settle  and  the  change  in  the  .level  of 
water  gave  the  desired  resin  volume.  The  density  was  then  calculated 
by  dividing  the  mass  of  the  oven-dry  resin  by  the  volume  of  the  wet 
resin.  The  density  was  measured  for  all  three  size  ranges.  The 
results  are  as  follows: 

Mesh  size  Density  gms*ml-‘L 


20,  50 
50,100 
100,200 


0.7988 

0.7975 

0.7937 


mean 


0.7960 
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V-B-2  Reactants: 

The  ethyl  acetate  used  was  99* *997  °]o  pure  with 
B.P.  =  76.9  -  77°  C 

Specific  gravity  =  0.894 

Molecular  Weight  =  88.107 

It  was  obtained  from  the  Fischer  Scientific  Co. 

Distilled  water  of  an  average  pH  of  6.8  was  used  in  the 
reaction.  This  was  mainly  to  prevent  the  degradation  of  the  resin. 

V-C  EXPERIMENTAL  PROCEDURE 

V-C-l  Measurement  of  Ester  Conversions; 

£alibration  of_flometers:  The  calibration  of  the  rotameters 

for  Ethyl  acetate  and  water  and  the  controlled  volume  pump  for  the 
effluent  stream  are  presented  in  FIGURES  (B-l)  to  (B-4)  of  Appendix  B. 

An_exp_erijii£ntal  run:  A  typical  experimental  run  to  measure  the 

ester  conversion  proceeded  as  follows. 

•  The  storage  tanks  for  distilled  water  and  ethyl  acetate 
were  filled,  whenever  necessary,  to  ensure  an  ample  supply 
of  the  reactants  for  3  to  4  runs.  At  the  start  of  a  run, 
the  feed  pumps  were  put  on  to  fill  the  overhead  feed  tanks. 

•  A  known  quantity  of  the  desired  size  of  the  oven- dry  resin 
was  introduced  in  the  reactor  followed  by  two  liters  of 
distilled  water.  The  reactor  was  covered  and  the  liquid 
level  inside  was  noted.  The  reactor  was  then  placed  in 

the  water  bath,  which  was  maintained  at  the  desired  temperature. 


, 
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Inlet  and  outlet  connections  were  made  and  the  reactor 
stirrer  was  started  and  adjusted  to  the  desired  speed. 
Within  a  short  time  the  reactor  temperature  reached  the 
bath  temperature.  Also  in  the  meantime  the  feed  pumps 
filled  the  overhead  feed  tanks  to  the  constant  .level  of 
the  overflow  lines. 

The  effluent  stream  pump  was  primed  and  its  stroke  length 
was  adjusted  for  the  desired  flow.  To  make  up  for  the 
depletion  of  the  liquid  level  in  the  reactor  caused  by 
the  priming  of  the  effluent  stream  pump,  fresh  distilled 
water  was  introduced  into  it  from  the  feed  tank. 

The  needle  valves  in  the  feed  lines  were  then  opened  to 
permit  the  flow  of  the  reactants.  Simultaneously,  the 
effluent  stream  pump  was  started.  The  rotameters  were 
quickly  adjusted  to  the  desired  flow  rates  of  the  reactants. 
The  reaction  had  now  started  and  was  allowed  to  reach  a 
steady  state.  During  the  transition  period,  steady 
conditions  of  flow  and  temperature  were  maintained  and 
the  stirrer  speed  was  frequently  checked.  A  constant 
liquid  level  in  the  reactor  indicated  that  inflow  was 
equal  to  outflow.  Alter  about  three  reactor  volumes  had 
passed  through,  samples  of  the  effluent  stream  were 
collected  and  analysed  for  the  ester  conversion.  The 
analysis  was  continued  till  three  consecutive  samples 
taken  at  intervals  of  one  half  the  residence  time  in  the 
reactor  agreed  within  +  Vjo. 
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Ana.l^;t_ical_methodsj_  The  analysis  of  the  effluent  stream  consisted 
in  measuring  the  acetic  acid  produced  in  the  reaction.  One  mole  of 
acetic  acid  was  formed  for  each  mole  of  the  ester  reacted.  No  attempt 
was  made  to  analyse  for  the  ethanol  and  the  unreacted  ethyl  acetate. 
Twenty-five  mis  of  the  effluent  stream  were  titrated  against  sodium 
hydroxide  of  known  strength  (xN)  using  phenolphthalein  as  an  indicator. 
The  volume  of  xN  NaOH  used  was  converted  to  an  equivalent  amount  of 
l/40  N  NaOH  for  consistency  in  reporting  the  data.  The  NaOH  solutions 
used  were  frequently  standardised  by  titrating  them  against  a  standard 
solution  of  oxalic  acid. 

Careful  checks  were  made  to  ascertain  the  extent  of  auto¬ 
catalysis,  if  any,  of  the  reaction.  The  acetic  acid  produced  in  the 
reaction  in  turn  may  catalyse  the  reaction.  Data  were  not  available  in 
the  literature  for  the  rate  of  hydrolysis  of  ethyl  acetate  catalysed  by 
acetic  acid.  However,  the  rate  was  experimentally  found  to  be  negligible 
for  the  residence  times  and  reaction  conditions  employed  in  the  study. 
This  was  also  observed  previously,  by  Harned  and  Pfanstiel  (26). 

It  was  observed  experimentally  that  distilled  water  (pH=6.8) 
kept  in  contact  with  the  resin  approached  a  constant  acidity  of  ph=4.1. 
The  final  pH  was  found  to  be  independent  of  the  resin  size,  temperature 
and  the  resin:  water  ratios  employed  in  the  experiments.  The  water  in 
contact  with  the  resin  quickly  approached  the  pH  of  4.1  and  remained 
constant  thereafter.  This  observation  was  further  confirmed  in  a 
continuous  system.  Distilled  water  was  passed  through  a  resin  bed  at  a 
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constant  flow  rate  (residence  time  in  the  bed  -  5  minutes)  for  about 
1.5  hours.  Throughout  this  period  the  pH  of  the  outgoing  water  remained 
steady  at  4.1.  Smith  and  Amundson  (54)  noted  the  same  phenomenon  and 
applied  necessary  corrections  in  their  analysis. 

This  increase  in  the  acidity  of  the  water  in  contact  with  the 
resin  can  be  attributed  to  the  presence  of  minute  quantities  of  metallic 
ions  (as  salts)  in  the  distilled  water.  It  is  believed  that  these 
cations  replace  the  H+  ions  in  the  resin  phase  and  the  H+  ions  in  turn 
pass  into  the  aqueous  phase  causing  a  decrease  in  the  pH.  This  also 
explains  why  the  pH  of  the  resin-treated  distilled  water  was  independent 
of  all  the  variables  mentioned  in  the  previous  paragraph. 

The  above  process  of  ion  exchange  is  accompanied  by  a  decrease 
in  the  activity  of  the  resin.  However,  calculations  showed  that  during 
an  average  service  of  a  batch  of  resin,  its  catalytic  activity  would 
decrease  by  about  1.1$.  This  deactivation  was  not  considered  to  be 
significant  and  the  resin  activity  was  assumed  to  be  constant  over  the 
duration  of  its  use.  The  day-to-day  variations  in  the  quality  of 
distilled  water  were  minor. 

The  acid  water  of  pH-4. 1^ was  used  to  make  a  5$  solution  of 
ethyl  acetate.  Ho  conversion  of  the  ester  was  noticed  even  after  four 
hours  of  stirring  at  100°  F. 
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In  the  analyses  of  the  experimental  data  it  was  therefore 
assumed  that  the  acetic  acid  produced  in  the  reaction  was  due 
solely  to  the  catalytic  effect  of  the  ion  exchange  resins.  To 
account  for  the  acidity  of  the  water  in  contact  with  the  cationic 
resins,  0.52  mis  of  l/40  N  NaOH  were  subtracted  from  the  readings 
for  the  resin  catalysed  reaction.  This  quantity  was  sufficient  to 
neutralize  25  mis  of  the  acid  water  of  pH  =  4.1. 

V-C-2  Residence  Time  Distribution  in  the  Reactor: 

With  a  view  to  evaluate  the  parameter  m,  the  residence 
time  distributions  of  liquid  in  the  reactor  were  measured  for 
different  stirrer  speeds  by  introducing  a  step  change  in  the  feed 
concentration.  Initially  the  steady-state  tracer  concentration 
in  the  reactor  was  CN(o)  and  then  at  9=0,  the  concentration  in  the 
inflow  was  reduced  to  zero.  The  response  was  measured  to  this 
negative  step  function  as  the  system  approached  a  new  steady  state. 
If  CN(q)  is  the  tracer  concentration  in  the  outflow  as  a  function 
of  time  the  response  of  the  proposed  model  is  (4.l): 

m  -  ■<•>  ■ 

t 

where,  9  =  — 

t 

The  tracer  chosen  was  hydrocholoric  acid;  the  use  of  a 
salt  was  undesirable  because  it  would  have  replaced  the  hydrogen 
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ions  from  the  resin.  The  residence  time  distributions  were 
measured  with  a  known  quantity  of  resin  in  the  system.  It  was 
thought  that  the  presence  of  solid  particles  may  have  seme  influence 
on  the  flow  pattern. 

Measurement  _of_the_tracer  concentration :  The  effluent  concentra¬ 
tion  of  HC.1  was  measured  as  the  impedence  of  the  solution  between 
two  electrodes  placed  in  a  impedance  cell.  The  cell  consisted 
of  two  26-guage  platinum  wires,  l/8  inch  long,  placed  l/2  inch  apart 
in  the  same  plane  and  encased  in  a  cylindrical  lucite  tube  of  l/4 
inch  i.d.  The  electrodes  were  plated  with  platinum  black.  The 
construction  of  the  cell  permitted  a  free  flow  of  the  liquid  through 
the  cell. 

The  cell  was  calibrated  with  HC.1  solutions  of  known  strengths. 
Each  standard  HC.1  solution  was  maintained  at  80°  F  by  using  the 
water  bath  because  the  impedance  of  the  solution  was  highly  sensitive 
to  the  temperature.  The  solution  was  withdrawn  at  a  steady  flow  by 
the  effluent  pump;  the  flow  rate  had  no  detectable  effect  on  the 
Impedance  over  a  moderate  range.  An  impedance  bridge  which  supplied 
an  a.c.  at  .1000  c/s  was  used  for  the  measurement.  This  high  frequency 
prevented  the  generation  of  gases  on  the  electrodes.  The  calibration 
curves  of  the  impedance  versus  strength  of  HC1  solutions  is  shown 
in  TABLE  (C-.1)  and  FIGURE  (C-l)  of  Appendix  C. 
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The  response  was  studied  by  recording  the  variation  of 
impedance  at  a  location  near  the  discharge  from  the  reactor  after 
the  negative  step  function  was  introduced  in  the  feed  stream.  Prior 
to  each  run  the  resin  particles  were  kept  in  the  HC.1  solution  of  the 
initial  concentration  CN(o)  to  insure  that  the  liquid  in  the  particles 
was  the  same  as  the  feed  liquid.  The  reactor  was  filled  to  a  volume 
of  slightly  more  than  two  .litres.  The  effluent  stream  pump  was  then 
started  to  withdraw  the  liquid  from  the  reactor  at  a  predetermined 
rate.  This  solution  as  it  passed  through  the  cell  gave  the  initial 
concentration  CN(o).  When  the  liquid  level  dropped  exactly  to  the 
reference  level  for  two  liters,  the  valve  for  distilled  water  feed 
to  the  reactor  was  quickly  opened  and  the  flow  adjusted  to  the  desired 
value.  Considering  the  relatively  high  residence  time  in  the  reactor 
-  .15  to  4 5  minutes  -  the  small  delay  in  this  procedure  was  assumed 
to  have  a  negligible  effect  on  the  response. 

The  transport  lag  in  the  system  was  taken  into  account.  A 
fluid  element  leaving  the  reactor  had  to  pass  through  nine  mis  of 
transport  section,  consisting  of  the  filter  and  ^  inch  stainless 
steel  tubing,  before  it  entered  the  impedance  cell.  It  was  assumed 
that  there  was  no  backmixing  in  this  section.  Therefore,  9/q 
minutes  after  the  step  change  was  imposed  (where  q  is  the  flow  rate 
of  the  effluent  stream  in  mls/min) ,  the  time  was  noted  as  9=0  and 
the  observations  started.  The  impedance  of  the  solution  was  then 
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noted  at  fixed  time  intervals.  This  was  continued  for  about 
three  residence  times  for  each  run. 

V-D  EXPERIMENTAL  RESULTS 

V-D-.1  Residence  Time  Distribution 

The  residence  time  distributions  were  measured  for  three 
different  conditions.  The  experimental  observations  are  presented 
in  TABLES  (C-2)  to  (C-4)  of  Appendix  C. 

The  experimentally  observed  quantities  were  9  and  the 
impedance  IMP(©).  The  corresponding  tracer  concentration  CN(©) 
were  read  from  the  calibrations  curve,  FIGURE  (C-l).  The  quantity 
=  l(©),  when  plotted  versus  9  on  semilog  graph  paper 

should  give,  according  to  equation  (V.2),  a  straight  line  of  slope=m 
and  intercepts.  These  plots  are  shown  in  FIGURES  (V-3)  to  (V-5)  . 

The  results  and  experimental  conditions  are  shown  in  TABLE  (V-l) . 

The  results  show  that  for  all  the  conditions  studied,  the 
parameter  m=l,  indicating  that  nearly  perfect  backmixing  occurs  in 
the  reactor  under  those  conditions. 

V-D-2  Significance  of  Liquid  Film.  Mass  Transfer 

Two  sets  of  runs  were  carried  out  to  evaluate  the  importance 
of  .liquid  film  mass  transfer  on  the  overall  conversion.  Turbulence 
in  the  flow  field  was  controlled  by  varying  the  stirrer  speed.  In 


CN(9) 

cn(o) 


FIGURE  V  ”  3  *  REACTOR  RESPONSE  TO  A  STEP  CHANGE  IN  THE  FEED 

CONCENTRATION 


FIGURE  V-4  * 


REACTOR  RESPONSE  TO  A  STEP  CHANGE  IN  THE  FEED 

CONCENTRATION 


FIGURE  V  -  5  : 
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TABLE  V-l 

EXPERIMENTAL  CONDITIONS  DURING  THE 
INVESTIGATION  OF  NON- IDEAL  REACTOR  OPERATION  AND  RESULTS 


Set  No. 

1-A 

1-B 

1-C 

N 

r  .p.m. 

75 

450 

55 

t 

min 

45 

45 

l4.8 

q. 

mls*min_1 

44.5 

44.5 

135 

¥ 

gms 

97 

97 

97 

T 

o  p 

80 

80 

80 

Fraction  of  the  feed 
entering  perfectly 

1 

1 

1 

backmixed  region  m 


TABLE  V-3 

TEMPERATURE  DEPENDENCE  OF  k  AND  D 


T 

1  x  103 

w 

k  x  102 

D  x  105 

°F 

°K 

T 

°K_1 

JL 

cm"2 

•  -1 
mm 

cm2/min 

80 

300.00 

3.3333 

36.52 

3.997 

2.9967 

100 

311.11 

3.2142 

4o.6o 

10.206 

6 .  .1915 

120 

322.22 

3.1034 

44.98 

23.726 

II.7272 
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each  set  the  turbulence  was  varied  over  a  wide  range  and  its  effect 
on  the  conversion  measured.  If  a  significant  film  mass  transfer 
resistance  existed,  a  change  in  the  conversion  with  increase  in  the 
stirrer  speed  would  be  noticed. 

A  known  amount  of  the  catalyst  was  introduced  in  the  reactor. 

The  resin  of  the  smallest  size  was  used  so  that  the  mass  transfer 

resistance,  if  any,  would  be  predominantly  evident.  The  flow  rate, 
q_,  and  temperature,  T,  were  maintained  constant  throughout  a  set, 
in  which  the  only  variable  was  the  stirrer  speed.  At  each  speed, 
the  conversion  of  ethyl  acetate  was  measured  after  a  steady  state 
was  reached.  The  experimental  conditions  for  the  two  sets  and  the 
results  are  presented  in  TABLES  (C-5)  and  (C-6)  of  Appendix  C. 

The  results  are  plotted  in  FIGURE  (V-6) .  A  reference  to  the  graphs 

shows  that  in  both  the  sets,  there  is  a  sharp  increase  in  the  conver¬ 

sion  from  a  stirrer  speed  of  100  to  175  r.p.m.  It  was  observed 
during  experiments  that  at  speeds  lower  than  100  r.p.m.  the  resin 
particles  were  completely  settled  on  the  reactor  bottom.  For  stirrer 
speeds  between  100  and  175  r.p.m.  the  particles  were  partially  in 
suspension,  the  degree  of  suspension  increasing  with  increase  in 
the  stirrer  speed,  and  for  more  than  200  r.p.m.  they  were  almost 
completely  suspended  in  the  bulk  of  fluid.  At  this  speed,  however, 
the  particles  were  not  uniformly  distributed  throughout  the  reaction 
mixture.  With  increase  in  agitation,  the  uniformity  in  the  catalyst 
distribution  increased  and  for  more  than  275  r.p.m.  the  particles 
were  completely  mixed  in  the  bulk. 
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FIGURE  V ”  8  :  EVALUATION  OF  THE  FILM  WSASS  TRANSFER 

RESISTANCE 


6o 


The  effect  of  stirrer  speed  on  the  mass  transfer  coefficient 


between  0-200  r.p.m.  is  therefore  complicated  by  the  unavailability 
of  all  the  catalyst  for  conversion.  The  conversion  after  the 
sharp  increase  in  that  range  of  stirrer  speed,  remains  almost 
constant  thereafter.  It  is  therefore  concluded  that  there  is  no 
significant  mass  transfer  resistance  for  the  stirrer  speeds  greater 
than  200  r.p.m.  Wo  conclusion  can  be  drawn  about  the  agitation 
range  of  0-200  r.p.m. 

In  the  second  set  of  observations  a  lesser  quantity  of  the 
catalyst  was  used.  This,  however,  made  no  difference  in  the  results. 
In  both  the  sets  the  temperature  in  the  reactor  was  80°  F.  Four 
additional  runs  were  made  at  100°  F  between  100  and  500  r.p.m.  No 
effect  of  agitation  rate  on  the  conversion  was  noticed  above  200 
r.p.m.  These  results  are  also  shown  in  TABLE  (C-7)  of  Appendix  C. 

V-D-3  Evaluation  of  Rate  Constants 

Having  established  that  under  the  conditions  studied  for 
the  system,  nonideal  flow  did  not  exist  (m=l)  and  the  film  mass 
transfer  resistance  was  negligible  (K  =  Go)  ,  equation  (lV..l4) 
reduced  to 


(IV..18) 


where, 


Q  _  Coth(wR)  -  l) 

■d2o 


R  ^ 


(iv. 15) 


w2  =  * 


D 


J  i 


' 


6.1 


This  linear  equation  permits  the  evaluation  of  the  rate  constants, 
k  and  D,  from  experimental  results. 

Four  sets  of  experiments  were  performed  at  a  constant 
temperature  of  80°  F,  with  different  but  known  amounts  (W)  and 
sizes  (R)  of  the  catalyst.  For  each  set,  the  ester  conversions 
were  measured  at  different  flow  rates,  q.  During  all  these  runs 
the  stirrer  speed  was  600  r.p.m.  The  experimental  data  are  given 
in  TABLES  (C-8)  to  (C-ll)  of  Appendix  C.  Each  set  of  runs  were 
plotted  as  _ 1  versus  q.  The  plots  are  shown  in  FIGURE  (V-7) • 

Ao 

The  functional  relation  between  the  variables  is  linear  and  the 
best  straight  .lines  for  the  experimental  data  of  each  set  were 
found  by  the  .least  square  method.  These  lines  were  forced  to  pass 
through  the  point  (0,  l)  as  it  should  according  to  the  equation 
(IV..18).  The  slopes  of  these  lines  is  equal  to  1  .  Thus  the 
values  of  Q,  were  obtained  with  the  corresponding  values  of  W 
and  R.  In  the  expression  (IV.15)  for  Q,  the  two  unknown  parameters 
are  k  and  D.  By  trial  and  error,  the  best  values  of  the  parameters 
were  calculated,  which  gave  the  minimum  standard  deviation  for  the 
experimental  data.  The  results  are  shown  in  TABLE  (V-2) .  They 
indicate  that  the  least- square-lines  are  approximated  very  closely 
by  the  theoretical  equations.  Details  of  the  calculation  procedures 
are  given  in  Appendix  C. 
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RATE  CONSTANTS 
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k  _and._  D  _as  f)ancti_on_of  t emp_er at ur _e :  The  effect  of  tempera¬ 

ture  on  k  and  D  was  investigated  by  also  determining  the  parameters 
at  100°  F  and  120°  F.  The  experimental  procedure  was  the  same  as 
above.  At  each  temperature  three  sets  of  runs  were  performed.  For 
each  set  W,  R  and  N  were  fixed  and  the  conversions  were  measured 
as  a  function  of  the  flow  rate,  q. 


TABLES  (C-12)  to  (C-17)  show  the  experimental  conditions 

and  data.  The  plots  of  versus  q  are  shown  in  FIGURES  (V-8) 

A 

Ao 

and  (V-9)  for  100°  F  and  120°  F,  respectively.  The  results  are 
given  in  TABLES  (V-2)  and  (V-3) . 

Both  k  and  D  were  found  to  be  exponential  functions  of 
temperature.  This  type  of  temperature  variation  for  the  diffusivity, 
D,  appears  to  be  a  characteristic  of  the  resin  catalyst  (52). 


The  parameters,  k  and  D,  are  plotted  versus  i.  in 
FIGURE  (V-10)  to  obtain  the  following  exponential  relations. 


l°g10k 

or  k 

and  log  D 

or  D 

where  T  is  in  °K . 


338^  +  9 . 88.15 

T 

5.665  x  1C)23  exp  ( -  min  ^ 

RT 

-  +  4.1577 

T 

3.747  x  .10^  exp(-  119.2Q)  cm2 -min 

RT 


-V 
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FIGURE  V-  8  : 


FIGURE  V  -  9  * 
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FIGURE 


V  -  10 


TEMPERATURE  DEPENDENCE  OF 


K  and  D 


c  ; 


CHAPTER  VI 


DISCUSSION 

Under  the  experimental  conditions  employed,  no  nonideality 
in  the  reactor  operation  was  observed.  This  observation,  in  retros¬ 
pect  does  not  seem  surprising  and  is  in  agreement  with  that  of  other 
workers  (59>6.l).  The  process  of  mixing  for  not  very  viscous  liquids 
has  been  shown  to  be  a  rapid  phenomenon  (7*38)*  In  an  ideal  CFSTR 
the  material  near  the  exit  is  moved  back  by  the  stirrer  to  the 
neighbourhood  of  the  entrance  in  a  time  much  less  than  t.  There¬ 
fore,  the  requirement  for  good  mixing  is  that  the  frequency  of 
rotation  of  the  stirrer  must  be  much  greater  than  l/t  .  This 
condition  was  amply  satisfied  under  the  present  experimental  conditions 
and  hence  no  apparent  difference  in  the  residence  time  distributions 
at  varying  speeds  was  observed. 

The  only  experimental  evidence  available  showing  the  actual 
nonideal  behaviour  of  a  CFSTR  is  the  work  done  by  Cholette  and 
Cloutier  (.10)  to  confirm  their  theoretical  model.  However,  they 
employed  a  specially  designed  vessel  which  had  such  features  that 
would  give  rise  to  dead  water  regions  and  by-pass  streams.  The 
reactor  had  an  i.d.  of  30  inches.  The  liquid  .left  the  vessel  by 
overflow  in  the  outlet  line  placed  to  give  a  liquid  depth  of  30  inches. 
The  feed  inlet  was  in  the  centre,  12  inches  from  the  surface  of  the 
liquid  and  the  agitator,  a  2.\  inch  turbine  model,  was  placed  in  between, 


■ 


■ 
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' 
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at  a  distance  of  6  inches  from  the  surface.  Apparently,  the 
positions  of  the  inlet,  outlet  and  the  stirrer,  contributed 
significantly  to  the  response  of  the  vessel.  They  obtained  m=0. 8 
at  150  r.p.m.  Their  results  however  would  be  applicable  only  to 
their  particular  set  up.  In  general,  it  can  be  concluded  that 
except,  possibly,  for  very  viscous  materials,  very  short  residence 
times  or  special  configuration  of  the  reactor  system,  the  agitator 
speed  has  little  influence  on  mixing.  Usually  other  factors  like 
heat  transfer  preclude  the  use  of  very  .low  agitation  rates. 

The  mass  transfer  resistance  in  the  external  film  was 
negligible  for  agitator  speeds  of  more  than  200  r.p.m.,  as  evidenced 
by  FIGURE  (V-6) .  The  reaction  rate  was  therefore  not  influnced  by 
the  film  diffusion  and  the  simultaneous  mass  transfer  and  chemical 
reaction  within  the  catalyst  particle  was  the  rate  limiting  step. 
Comments  on  the  effect  of  film  diffusion  at  agitator  speeds  less 
than  200  r.p.m.  is  not  possible  because,  in  that  region,  the  inter¬ 
facial  area  was  not  independent  of  the  agitation  rates.  Then, 
obviously,  the  total  amount  of  the  catalyst  was  not  measurably 
productive.  For  a  meaningful  analysis  of  such  a  phenomenon  the  rate 
must  be  expressed  per  unit  of  catalyst  on  which  the  measured  conversion 
occurs.  Mount  of  the  effective  catalyst,  however,  was  not  possible 


to  determine. 


‘ 
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Very  few  investigators  working  with  similar  systems  have 
observed  external  mass  transfer  as  a  rate  .limiting  process.  Siato 
et  al  (See  (49)  )  observed  deviation  from  the  first  order  kinetics 
of  the  inversion  of  sucrose  indicating  an  effect  of  external  diffusion. 
However,  the  later  investigation  of  Reed  and  Dan dr off  (49)  for  the 
same  reaction  system  indicated  no  effect  of  film  diffusion  over  a 
wide  range  of  Reynold  Nos.  Engel  and  Hougen  (l8)  observed  very  strong 
external  diffusion  effects  for  particles  of  0.00256  cms  diameter  at 
80°  F  and  580  r.p.m.,  while  for  0.01153  cms  diameter  at  126°  F  and 
175  r.p.m.  the  film  diffusion  was  negligible.  In  the  present  study 
the  external  mass  transfer  resistance  was  observed  to  be  negligible 
even  for  0.00283  cms  diameter  particles  at  100°  F.  The  descrepancy 
cannot  be  simply  explained.  However,  it  should  be  noted  that 
Engel  and  Hougen  worked  with  amyl  acetate  while  the  present  work 
was  concerned  with  ethyl  acetate. 

The  last  phase  of  the  study  was  concerned  with  evaluating 
the  parameters  k  and  D  and  their  dependence  on  temperature.  Initial 
observations  indicated  the  overall  reaction  rate  per  unit  mass  of 
catalyst  as  measured  experimentally  was  almost  independent  of  the 
particle  size  of  the  catalyst  and  hence,  was  limited  by  the  relatively 
slow  chemical  reaction;  the  internal  diffusion  had  no  significant 
influence  on  the  reaction  rate.  As  a  result  of  this  it  was  found, 
while  calculating  the  parameters  k  and  D  from  the  experimental  data, 
that  for  small  particles  the  diffusivity  can  be  varied  over  relatively 
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wide  limits  without  markedly  affecting  the  value  of  Q.  However,  the 
effect  of  diffusivity  increased  with  the  increase  in  the  particle  size. 

The  effectiveness  factors  calculated  from  equation  (II.5) 
for  different  particle  sizes  at  different  temperatures  are  shown  in 
TABLE  (VI-l)  and  FIGURE  (VI-l) . 

For  the  ion  exchange  resin  used  in  this  study,  the  diffusivity 
was  also  found  to  be  an  exponential  function  of  temperature,  though 
less  strong  than  the  reaction  rate  constant.  The  parameters  w  and 
are  therefore  only  weakly  dependent  on  temperature.  With  other 
solid  catalysts  for  which  the  internal  diffusivity  is  not  an  exponential 
function  of  temperature,  the  effectiveness  factor,  ,  would  be 
expected  to  decrease  more  rapidly  and  the  overall  reaction  rate  to 
increase  more  slowly  with  increase  in  temperature. 

A  reference  to  TABLE  (Vl-l)  shows  that  for  the  resin  particles 
below  0.0254  cms  in  diameter  and  120°  F  the  effectiveness  factors  are 
very  near  to  unity  and  the  rigorous  rate  expression. 

Rate  =  T^kcK 
can  be  approximated  by 

Rate  =  kC^ 

$ 


with  deviations  of  only  up  to  2.1 % 


‘  ; 
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TABLE 

VI-1 

EFFECTIVENESS 

FACTORS 

T 

OJP 

R 

cms 

w 

cm”1 

WR 

80 

0.03045 

0.01268 

0.00283 

36.52 

0.1112 

0.0450 

0.0103 

100 

0.03045 

0.01268 

0.00283 

4o.6o 

0.1238 

0.0515 

0.0115 

120 

0.03045 

0.01268 

0.00283 

44.98 

O.I37O 

0.0570 

0.0127 

n 


0.9262 

0.9860 

0.9993 

0.9110 

0.9827 

0.9991 

0.8938 

0.9789 

0.9989 
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FIGURE  VI -I  *  EFFECTIVENESS  FACTORS 
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The  equations  (A. 9)  and  (A.IO)  give 


C 

C 


s 


C 

C 


o 


R  Sinh 
r  Sinh 


from  which  the  internal  concentration  profile  for  the  ester  in  terms 
of  the  surface  concentration  (which  in  this  case  is  the  effluent 
concentration)  can  be  calculated.  The  results  of  these  calculations 
are  shown  in  TABLE  (VI-2)  and  FIGURE  (VI-2) .  The  flat  nature  of  the 
curves  indicate  relatively  fast  internal  diffusion. 


The  high  effectiveness  factors  provide  an  indirect 


confirmation  of  the  observation  that  the  external  film  diffusion  is 
not  a  rate  .limiting  step.  As  Petersen  (48)  says, "Resistance  to  mass 
transfer  within  a  catalyst  is  greater  than  the  resistance  between 
the  external  surface  and  the  bulk  fluid  phase.  This  means  that  if 
diffusion  to  the  external  surface  limits  the  overall  reaction  rate, 
then  internal  diffusion  must  also  limit  the  reaction  rate  and  the 
designer  is  obliged  to  make  an  analysis  which  includes  both  external 
and  internal  diffusion.  One  can  therefore  dismiss  the  practical 
possibility  of  a  system  in  which  external  diffusion  is  controlling 
in  coexistance  with  an  effectiveness  factor  of  unity. " 


Since  the  above  work  was  finished,  Bochner  et  al  (8) ,  who 


studied  the  esterification  of  salicylic  acid  and  methanol  using 
Dowex  50\f -x8  resins  as  the  catalyst,  observed  that  the  external  and 


internal  diffusion  of  the  reaction  species  had  no  significant 
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FIGURE  VI- 2  *  CONCENTRATION  PROFILE  OF  ETHYL 


ACETATE  IN  THE  RESIN  PARTICLES 
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influence  on  the  overall  reaction  rate.  They  demonstrated  that  the 
competitive  chemisorption  of  water  and  salicylic  acid  was  the  rate 
limiting  step  and  a  Langmuir  type  equation  fitted  the  date  for 
0  -  84.5  °!o  conversion. 

It  was  not  the  aim  of  the  present  study  to  evaluate  such 
models  and  formulate  the  actual  reaction  mechanism.  The  ester  con¬ 
centrations  and  the  conversions  in  the  reaction  were  small.  It  is 
possible  that  the  ester  concentration  within  the  catalyst  volume 
would  be  higher  than  that  in  the  surrounding  liquid.  However,  a 
linear  distribution  law  may  correctly  relate  concentrations  in  the 
bulk  liquid  and  pore  liquid  (8).  In  such  a  case,  then,  the  calculated 
chemical  reaction  rate  constant,  k,  would  be  the  product  of  the 
distribution  coefficient  and  the  true  reaction  rate  constant. 

From  FIGURE  (V-,10)  the  value  of  the  rate  constant  is 
k(resin)  =  0.04.1  min  ^  at  25°  C 

For  an  equivalent  H+  ion  concentration  the  value  of  k 
for  the  hydrolysis  of  ethyl  acetate  in  a  homogeneous  medium  using 
HC.1  acid  is  (26): 

k(HCl)  =  0.0409  min-1  at  25°  C 

5  1.0  at  25°  C 


k(resin) 

k  (HC.1) 


I 
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Bernhard  and  Hammett  (5)  obtained  k(resin)/k(HCl)  =  O.92 

at  25°  C,  which  is  in  reasonable  agreement  with  the  above  value. 
Saletan  and  White  (52)  obtained  the  ratio  equal  to  1.84  for  the 
esterification  of  acetic  acid  and  ethanol  at  25°  C.  The  value  of 
E,  the  energy  of  activation  for  the  reaction,  however  agrees  more 
closely  with  that  of  Davies  and  Thomas  (l6)  as  shown  below: 


E  = 

15,500 

cal*  gmmo.le"-1' 

(present  work) 

E  = 

17,500 

cal-gmmole"1 

(Bernhard  and  Hammett  (5)  ) 

E  = 

15,540 

cal  •gmmo.le”"1- 

(Davies  and  Thomas  (l6)  ) 

The  following  comparison  of  the  values  of  the  effective  diffusivity, 

D,  shows  good  agreement  in  the  order  of  magnitude. 

Ethyl  Acetate  D  =  3.1  x  10"  ^  cm^-min-^  at  25°  C  (present  work) 

Ethyl  Formate  D  =  5-0  x  10” 5  cm2. min”-1-  at  25°  C  (Smith  and  Anundson(54' 

Ethanol  D  =  1.5  x  10" ^  cm^-min'-1-  at  25°  0  (Saletan  and  White (52)) 

The  "activation  energy"  for  the  diffusivity  D  of  ethyl 
acetate  in  Dowex  5OWX8  was  found  to  be  11,920  cal* gmmo.le"”1’  (FIGURE 
(V-IO)  )  as  compared  to  that  of  ethanol  in  Dowex  5O-X8  equal  to  9 >850 
cal-gmmole"1  (.12).  In  .light  of  the  fact  that  different  types  of 
resins  were  employed  by  different  workers,  further  comparison  is 
not  pertinent.  However,  the  general  agreement  with  the  previous 
works  in  the  order  of  magnitudes  of  k  and  D  and  the  relatively  small 
scatter  in  the  fit  of  the  experimental  data  to  the  proposed  equations 
provides  a  good,  though  only  partial,  verification  of  the  model. 


, 
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It  should  be  noted  that  the  effective  diffusivity,  D,  is 
different  than  the  true  diffusional  coefficient.  In  reality  the 
diffusiona.1  process  in  solids  is  an  impressively  complex  matter. 

The  effective  diffusivity  takes  into  account  the  porosity  of  solid. 

In  addition,  all  the  pore  volume  may  not  be  effective  and  different 
portions  of  the  surface  are  characterized  by  different  accessibilities. 
The  transport  within  the  solid  involves  flow  splitting  and  covergence 
and  multidirectional  flow  paths  of  which  the  individual  contributions 
to  the  flow  resistance  are  unknown.  One  may  lump  all  these  factors 
into  some  parameters;  they  are  calculational  conveniences  and  do  not 
throw  any  direct  light  on  the  flow  mechanism.  In  view,  therefore, 
of  the  limitations  attendant  upon  employing  any  sophisticated  model, 
further  refinement  hardly  seems  profitable  from  the  engineering 
point  of  view.  The  use  of  effective  diffusivity  enables  one  to 
describe  the  diffusion  process  independent  of  the  shape,  size  and 
distribution  of  the  pores. 

A  complete  verification  of  the  mathematical  model  is  still 
to  be  achieved.  The  absence  of  any  significant  diffusional  resis¬ 
tance  in  the  reaction  and  any  nonideality  in  the  liquid  flow  through 
the  reactor  left  only  a  simple  model  to  be  verified.  If  the  proposed 
model  were  realistic,  one  would  expect  to  obtain  m=.l  for  all  the 
experimental  runs  performed  at  high  stirrer  speeds.  In  other  words, 
the  best  straight  line  passing  through  the  experimental  points  on 
the  plots  of  Q,  versus  q  (FIGURES  (V-7)  to  (V~9)  )  should  intercept 
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the  vertical  axis  at  unity.  The  values  of  the  intercepts  thus 
obtained  from  the  experimental  data  were  found  to  be  varying  between 
0.7  to  1.2.  It  is  possible  that  the  proposed  model  is  inadequate. 
However,  it  is  thought  that  more  experimental  runs,  particularly 
at  low  flow  rates,  q,  should  be  performed  for  a  meaningful  evalua¬ 
tion  of  the  model.  This  was  not  feasible  with  the  present  set  up 
due  to  the  equipment  limitations. 

A  necessary  condition  for  a  meaningful  evaluation  of  the 
complete  model  is  that  the  experiments  should  be  carried  out  at 
agitator  speeds  at  which  the  catalyst  particles  are  completely 
suspended  in  the  reactor  volume.  The  intensity  of  turbulence  at 
these  speeds  is,  however,  usually  sufficient  to  eliminate  nonideality 
in  the  flow  patterns  and  the  external  film  resistance.  While  not 
much  can  be  done  about  the  nonideality  in  flow,  it  is  possible  that 
by  judicious  choice  of  certain  factors  the  external  mass  transfer 
resistance  can  be  significantly  increased.  For  example  even  when 
the  particles  are  suspended  in  the  bulk  of  liquid,  the  relative  motion 
between  the  two  phases  may  be  decreased  by  eliminating  baffles  and 
using  an  appropriate  agitator.  A  higher  ester  may  be  employed.  Also 
since  the  reaction  rate  increases  exponentially  with  temperature, 
while  mass  transfer  increases  little  more  than  linearly  one  would 
anticipate  that  at  sufficiently  high  temperatures  the  rate  might 
become  mass  transfer  controlled  for  all  catalysts. 


8.1 


It  is  theoretically  possible  to  incorporate  other  modes 
of  flow  in  a  mathematical  model,  but  the  quantitative  expressions 
become  too  complex  to  be  of  any  practical  use.  For  example, 
accounting  a  fraction  of  the  reactor  as  a  plug-flow  region  would 
introduce  exponential  terms  in  the  model  equation  and  the  parameters 
in  this  equation  cannot  be  easily  determined  by  experiments. 

Presence  of  a  dead  water  region  in  the  vessel  where  there 
are  no  solids  does  not  affect  the  conversion.  Finally,  the  problem 
becomes  very  difficult  when  the  solid  particles  are  settled,  at  the 
reactor  bottom.  Then  a  fraction  of  the  solids  which  have  settled 
is  also  effective  in  the  reaction.  The  total  amount  of  catalyst 
effective  in  this  case  remains  virtually  indeterminate  as  noted 
previously. 


' 
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NOMENCLATURE 

Ephemeral 

notation  whose  bearing  does  not  extend  beyond  the 

immediate 

context  is  not  included. 

Ao 

Concentration  of  acid  in  product  stream,  gmmoles  *m.l~^ 

C 

Concentration  of  ester  in  the  resin  particle, 

gmmoles  •ml""'" 

ti 

•H 

O 

Concentration  of  ester  in  feed,  gmmoles 'ml"1 

o 

o 

II 

Concentration  of  ester  in  product  stream,  gmmoles 'ml--1- 

Cs 

Concentration  of  ester  on  the  external  surface  of  the 

resin  particle,  gmmoles *ml“^ 

c* 

Concentration  of  ester  in  exit  stream  of  the  perfectly 

CN(©) 

backmixed  part  of  a  CFSTR 

Tracer  concentration  at  time  0,  gmmoles  •ml""*' 

D 

Effective  diffusivity  of  ester  in  the  resin  phase, 

2  •  -1 
cm  -mm 

E 

Energy  of  activation,  cal'gmmole-1' 

F(0) 

Residence  time  distribution  function 

IMP 

Impedance  of  HC.1  solution  in  the  .lucite  cell,  ohms 

K 

Mass  transfer  coefficient  for  ester  diffusion  through 

liquid  film,  cms-min"1;  K  ;  at  terminal  speed 

N 

Stirrer  speed,  r.p.m. 

Q 

3^(wR- coth(wR)  -  l) ,  mls-min"1 

R2^ 

R 

Radius  of  resin  particle,  cms 

Rv 

Reaction  rate  per  unit  volume  of  the  catalyst,  gmmoles •m.l-1 

per  minute 


.  u  V 
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T 

V 


W 

X(lOO) 

d 


k 

m 

d 


r 


t 

t 

w 

z 

0 

n 

? 


Temperature 

Volume  of  N/40  NaOH  required  to  neutralize  acetic  acid 
in  25  mis  of  product  stream,  mis 
Weight  of  resin  catalyst,  gms 
Percentage  conversion  of  ester 

Diameter  of  resin  particle;  daVg  =  the  surface  area 
average  diameter  of  the  resin  particles,  cms 
Chemical  reaction  rate  constant,  min'-1- 
Fraction  of  feed  by-passing  the  catalyst  in  a  CFSTR 
Feed  rate;  1  =  ester  feed  rate,  d  =  water  feed 

rate,  mls-min“-L 

Space  variable  in  a  resin  particle;  distance  from  the 
centre,  cms 
Time,  min 

Average  residence  time  of  a  reaction  mixture  in  CFSTR,  min 


(wR-coth(wR)  -  l) 

Reduced  time,  t/t 

Effectiveness  factor  of  a  catalyst 

Density  of  resin  catalyst,  gms -ml-1 
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APPENDIX  A 


DERIVATION  OF  THE  MATHEMATICAL  MODEL 


The  linear  differential  equation 


d2C 


+ 


2/d C\  kC  = 


dr  r  'dr 


0 


D 


(A.l) 


has  a  general  solution  of  the  form  (See  CHAPTER  IV) 


0  _  A  cosh(wr) 


B  sinh(wr) 
r 


(A.  2) 


where, 


wL 


k 

D 


The  boundary  conditions  imposed  by  the  system  are: 
(l)  At  the  centre  of  the  particle 


(■ 


dC' 

drj 


=  0 


r=o 


(A. 3) 


The  equation  (A. 2)  now  reduces  to 


B  sinh(wr) 
r 


(A.  b) 


(2)  At  the  catalyst  surface  (See  equation  (CV.  13)) 


3WD  +  D 
R^mq  K 


dC\ 

Wr=R 


(A.  5) 
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Differentiating  equation  (A. 4) 


B  (wR- coth(wR)  -l) 

2  sinh(wR) 
r  v  ' 


(A. 6) 


Also, 


Cs 


(C)r=R 


B  sinh(wR) 
R 


(A.  7) 


Substituting  equations  (A. 6)  and  (A. 7)  in  equation  (A. 5)  and 
simplifying, 


B 


Ci 

(l  +  Pz) 


R 

sinh(wR) 


where 


P 


3WD  +  D_ 
R^mq  RK 


and 


(A.8) 


z  =  (wR* coth(wR)  -  l) 


Equations  (A.8)  and  (A.  4)  give 


C  1  R  sinhfwr^ 

Cj_  _  (l  +  Pz)  r  sinh(wR) 


(A.9) 
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At  the  particle  surface 

r  =  R 

-  —  1  (A.  10) 

Ci  r=R  Ci  1  +  Pz 

Now,  according  to  the  model  proposed 

C*  =  Bulk  concentration  of  the  ester  in  the 
perfectly  backmixed  part  of  the  reactor. 
Therefore,  a  rate  balance  for  the  diffusion  of  the  ester  through 
the  liquid  film  surrounding  a  particle  gives 


»  1 


r=R 


=  K(C*  -  Cs) 


Rearranging  the  above  expression 


* 


1  +  ”  d  (5 

)i  K  Ci  W, 


r=R 


(A. 11) 


From  equations  (A. 8)  and  (A. 6) 


dCN 

s.dr/  r=R 


Ci 


R  (1  +  Pz) 


(A. 12) 


Substituting  (A. 12)  and  (A.IO)  into  (A.ll) 


* 


1 


(1  +  Pz) 


D  z 
RK  (1  +  Pz) 


+ 


(A. 13) 
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However,  from  equation  (IV. 4) 


Cn  mC* 

—  =  .1  -  m  +  _ 

Ci  C± 


Cn  m 

—  =  1  -  m  +  _ 

Ci  (.1  +  Pz) 


1  + 


Dz 

RK 


Simplifying  the  above  expression 


z  3WD 

(1  +  Pz)  R^q 


Ao  _  ,  Op 

Ci  "  "  Ci 


3WD 


(1  +  Pz)  R^q 


Rearranging  equation  (A.l6) 


Ci 


A, 


o 


R%q.  +  Rfq  +  i 

3WDz  3KW  m 


or 


A  -  i 

A  „ 


1 

Q 


+ 


R^ 

3KW 


1 

+  _ 
m 


where 


Q  =  (wR'coth(wR)  -  l) 

R2? 


(A.,14) 


(A.15) 


(a.i6) 


(A.  17) 


(A.l8) 
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(A. 17)  and  (A.l8)  are  the  desired  expressions  for  the 
proposed  model. 

Effectiveness  Factor; 

When  K  =  00  and  m  =  1,  equation  (A.  .17) 


reduces  to 


and 


However , 


A  .  1  +  * 

Ao  Q, 


°s  =  c*  =  C0 


Ci  =  Co  +  Ao 


Lo  q 

i  +  —  =  i  +  - 
Aq  q, 


or 


qAQ  =  Q,C 


o 


The  rate  of  disappearance  of  ester 


(A. 19) 


(A.  20) 


(A.  21) 


=  qA0  gmmoles/min 

Therefore,  rate  of  the  reaction  per  unit  effective  volume  of 


the  catalyst  is  given  by 


. .  ;"jOL  O  '■  , 


92 


R 


V  = 


$<iA 


o 


W 


-1  -1 

gmmoles*min  •  ml  of  catalyst 


£QC 


o 


W 


(Using  equation  (A.2l)  )  (A. 22) 


Wow,  by  definition 


Rv  =  ^kc< 


=  Y|kCQ  (Using  equation  (A.  20)  )  (A.  23) 


where, 


T]  =  effectiveness  factor 
Comparing  equations  (A. 22)  and  (A. 23) 

=  i  S-  (wR-coth(wR)  -  1) 
W  kR^ 


—2 —  (wR*coth(wR)  -  l) 
w2R2 


(A. 24) 


or , 


No 

1 

4 

7 

10 

13 

l6 

19 

22 

25 

28 

31 

34 

37 

4o 

43 

46 

49 

52 

55 

58 

6i 

64 

67 

70 

73 

76 

79 

82 

85 

88 

91 

94 

97 

.100 
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APPENDIX  B 


MISCELLENEOUS  DETAILS 


TABLE  B-.1 


MEASUREMENT  OF  PARTICLE  DIMETER 


(a)  Size  Range:  20-50  mesh 


Diameter 

No. 

Diameter 

No. 

Diameter 

_  O 

d  x  10  c  cms 

d  x  10‘ 2  cms 

d  x  10" 2  cms 

1.6956 

2 

3.4916 

3 

6.0200 

4.4247 

5 

8.0969 

6 

4.5451 

5.1772 

8 

6 . 0200 

9 

4.4247 

6 . 3210 

11 

3.7625 

12 

5.7190 

4.4548 

l4 

6 . i4o4 

15 

4.7558 

5.8190 

17 

4.3043 

18 

6 . 1805 

4.6053 

20 

3.0401 

21 

6.9631 

6 . 0200 

23 

5.5384 

24 

8.2183 

6 . 3618 

26 

3.7023 

27 

2.4090 

6 . 4113 

29 

5.1772 

30 

3-7324 

5.3227 

32 

3-4013 

33 

2.3789 

4.6956 

35 

4.6655 

36 

4.4648 

5 . 9899 

38 

3.624.1 

39 

6 . 3210 

1.2642 

4l 

2.8595 

42 

4. 4648 

4.3344 

44 

2.7993 

45 

6 . 8026 

2.5294 

47 

5.1772 

48 

4.0635 

4.3650 

50 

6.6220 

51 

6 . i4o4 

6.1705 

53 

6.8628 

54 

6.3812 

4.0033 

56 

5.0879 

57 

7.0735 

7.5551 

59 

4.7558 

60 

6 . 5618 

4.0033 

62 

5.3578 

63 

6.7424 

5.5685 

65 

6 . 0802 

66 

7.2240 

4.2i4o 

68 

7.4949 

69 

4.8160 

4.5155 

71 

3.1906 

72 

7.1036 

7.1046 

74 

5.9899 

75 

1.8060 

5.1471 

77 

7 . 4o46 

78 

5.8695 

3.3110 

80 

5.9598 

81 

7.5250 

4.4849 

83 

5.5685 

84 

7.7056 

7.2842 

86 

5.8695 

87 

2.2274 

7.4648 

89 

3.0100 

90 

4.6354 

2.7692 

92 

4.7869 

93 

5.9297 

3.3712 

95 

4.7859 

96 

2.6187 

4.0033 

4.7257 

98 

7.5551 

99 

4.3645 

94 


100 

XI  dj2  =  0.2874 
1 


100 

'21  a.3 
1  1 


=  I.7503  x  10 


-2 


avg 


 Idi 


.3 


and,  Relative 

Standard  Deviation 


=  6.O9O  x  10" ^  cms 


R  = 


0.03045 
100 


-  davg) 


cms 

2 


(100  -  1)  d2ayg 


x  100 


=  28.7  i 


TABLE  B-l  (CONTINUED) 
MEASUREMENT  OF  PARTICLE  DIAMETER 


(l)  Size  Range:  50-100  mesh 


No. 

Diameter 

No. 

Diameter 

No. 

Diameter 

d  x  10" 2  cms 

-  P 

d  x  10  cms 

O 

d  x  10" ^  cms 

1 

1.4308 

2 

2.  5165 

3 

2.8472 

4 

2.5524 

5 

2.7178 

6 

1.3733 

7 

2.8832 

8 

2.5309 

9 

2.4086 

10 

1.8478 

11 

2.7260 

12 

1. 5746 

13 

2.1714 

l4 

2.6200 

15 

2.8688 

16 

1.8263 

17 

2. 5894 

18 

1.5243 

19 

2.5165 

20 

2.0995 

21 

2.5093 

22 

1.6893 

23 

2.5884 

24 

2.8185 

25 

1.9000 

26 

2 . 4o86 

27 

2.1858 

28 

2.1383 

29 

2.4662 

30 

1.9195 

31 

3.1825 

32 

2.7394 

33 

2.5956 

(CONTINUED) 
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TABLE  B-,1  (CONTINUED) 
MEASUREMENT  OF  PARTICLE  DIAMETER 


(b)  Size  Range:  50-100  mesh 


No. 

Diameter 

No. 

Diameter 

No. 

Diameter 

_  O 

d  x  10  c  cms 

d  x  10“  ^  cms 

d  x  10“ ^  cms 

34 

2.7681 

35 

2.5309 

36 

2.9361 

37 

2.6028 

38 

2.4374 

39 

2.7610 

40 

2.7466 

4i 

2.2864 

42 

1.7472 

43 

2.9048 

44 

2.6747 

45 

1.3373 

46 

3.5734 

47 

2.1642 

48 

2.6423 

49 

2.3367 

50 

2.4805 

51 

1.6321 

52 

1.8615 

53 

1.5530 

54 

3.0486 

55 

2.4158 

56 

2.7538 

57 

2. 5740 

58 

2.7034 

59 

2.8275 

60 

2.5524 

6l 

2.6603 

62 

2.4156 

63 

2.6387 

64 

2.6315 

65 

2.5596 

66 

2.3439 

67 

1.9772 

68 

1.8400 

69 

1.7975 

70 

2.7897 

71 

2.3799 

72 

1.6537 

73 

2.0923 

74 

2.4230 

75 

2.7466 

76 

2.  5940 

77 

2.6315 

78 

2.6459 

79 

2.5596 

80 

1.6968 

81 

2.3224 

82 

2.6028 

83 

1. 5027 

84 

1.8047 

85 

2.5956 

86 

2.7394 

87 

2.4374 

88 

2.4158 

89 

1.5099 

90 

3.0054 

91 

2.4158 

92 

2.1642 

93 

1.6171 

94 

1.3733 

95 

2.4805 

96 

2.3871 

97 

2.7034 

98 

1.8760 

99 

2.7897 

100 

2.0923 

100 

ZL  =  0.5766 


100 

XI  =  1.4623  x  10“3 
l 


avg 


z 


d,3 


•  * 


• 
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d. 


avg 


R  = 


2.536  x  10“ 2 

0. 01268 


cms 


cms 


and,  Relative 

Standard  Deviation 


100  2 

H(di  "  d-avg) 

1 


(100  -  1)  d 


avg 


=  19.5  i 

TABLE  B-l  (CONTINUED) 
MEASUREMENT  OF  PARTICLE  DIAMETER 


JL. 

2 


x  100 


(c)  Size  Range:  100-200  mesh 


No. 

Diameter 

No. 

Diameter 

No. 

Diameter 

d  x  10“ 2  cms 

d  x  .10" 2  cms 

O 

d  x  10  c  cms 

1 

6.0385 

2 

8.0014 

3 

5.4484 

4 

6 . 4240 

5 

2.1074 

6 

7.1703 

7 

1.6191 

8 

2.7756 

9 

6.7334 

10 

1.4906 

11 

7.4787 

12 

6.3993 

13 

1.4906 

l4 

5.4484 

15 

2.5957 

l6 

0.8995 

17 

2. 5186 

18 

4.1634 

19 

4.39^7 

20 

3.8550 

21 

1.3364 

22 

1.7733 

23 

7-3245 

24 

5.7054 

25 

1.0537 

26 

3.2125 

27 

0.8738 

28 

7 . 04l8 

29 

5.3713 

30 

2. 5700 

31 

5.4227 

32 

4.2405 

33 

2.7499 

34 

1 . 8247 

35 

4.8059 

36 

5.7054 

37 

3.3153 

38 

4.3934 

39 

7.8385 

40 

5.0372 

4.1 

3.3667 

42 

4.2262 

43 

4.8830 

44 

6 . 6007 

45 

5.9367 

46 

5.2428 

47 

1.4,135 

48 

2.2616 

49 

2.7499 

50 

2.1845 

51 

1.5163 

52 

4.8573 

53 

1 . 8761 

54 

1.6705 

55 

4.0349 

56 

3.8807 

57 

4.7031 

58 

2.6985 

59 

3.1354 

60 

5.5769 

61 

2.3387 

62 

1.9275 

63 

5.7285 

64 

4.3967 

65 

2.4415 

66 

7.9670 

67 

4.0349 

68 

7.016.1 

69 

2.2359 

( CONTINUED) 
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TABLE  B-l  (CONTINUED) 
MEASUREMENT  OF  PARTICLE  DIAMETER 


(c)  Size  Range:  100-200  mesh 


No. 

Diameter 

No. 

Diameter 

No. 

d  x  10- 2  cms 

d  x  10“ ^  cms 

70 

4.0763 

71 

2.7756 

72 

73 

1.3364 

74 

3.6237 

75 

76 

4.6603 

77 

2.2359 

78 

79 

6 . 6049 

80 

7.2731 

81 

82 

7. 5000 

83 

6.4764 

84 

85 

2.5957 

86 

2.3644 

87 

88 

5.3199 

89 

1.6705 

90 

91 

1.4649 

92 

5. 4484 

93 

94 

2.3901 

95 

5.3456 

96 

97 

1.5934 

98 

1. 5163 

99 

100 

3.9814 

100 

Z!  di  =  1.973  x  10-3 
1 


100  o  «- 

^  <^3  =  1.1167  X  10“  7 

1 


davg 


R 


and,  Relative 

Standard  Deviation 


=  5.68  x  10“ 3  cms 

=  0.002843  cms 

Too  ~2~ 

( d-i  ^avg) 

(lOO  -  l)  d2avg 


i_ 

2 


x  100 


Diameter 

x  10" 2  cms 

4. 5489 
5.0629 
6.0395 
1.4649 
1.5677 
5.2488 
3.3667 
5 . 4484 
1 . 4906 
2.5443 


37.2  i 
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TABLE  B-2 
EQUIPMENT  DETAILS 


Equipment 


(1) 

Storage  bottle  for 
distilled  water 

Material 

Size 

(2) 

Storage  bottle  for 
ethyl  acetate 

Material 

Size 

(3) 

Constant  head  tank 
for  water 

Material 

Size 

(4) 

Constant  level  tank 
for  ethyl  acetate 

Material 

Size 

(5) 

Feed  pump  for  water 

Material 

Type 

Size 

(6) 

Feed  pump  for  ester 

Material 

Type 

Capacity 

Plunger 

Discharge 

Pressure 

(7) 

Rotameters 
( a)  for  water 

Material 

Range 

Type 

(b)  for  ethyl 
acetate 

Material 

Range 

Type 

(8) 

Reactor 
(a)  Vessel 

Material 

Size 

Details 

Pyrex  glass 
12  gallons 

Pyrex  glass 

1  gallon 

glass 
1000  mis 

glass 
500  mis 

monel 

Centrifugal,  Fisher  Scientific  Cc 

2  inch  diameter  turbine  capable 
of  supplying  1  gpm  at  12  ft.  of 
water  head. 

Stainless  Steel 

Positive  displacement,  duplex, 
Milton  Roy  Co. 

5  mls/min/side 

l/lO  inch  dia.  x  1  inch  length, 

96  s/min. 

2000  psia 

glass 

O-.13O  mls*min"l 
No.  603,  Matheson  Co. 

glass 

0-5  mis  •min” 

No.  601,  Matheson  Co. 

Pyrex  glass 

9y  inches  height,  6  inches  o.d. 

5y  inches  i.d. ,  3/^  inch  wide 
flange,  ground  top. 


(CONTINUED) 
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TABLE  B-2 


(8) 


(9) 


EQUIPMENT  DETAILS  (CONTINUED) 


Equipment 

Details 

Reactor 

Material  : 

Stainless  Steel 

(b)  Baffles 

Size  : 

No.  : 

\  inch  width,  9ir  inches 
length,  l/8  inch  thickness. 

4 

(c)  Filter 

Material  : 
Size  : 

Sintered  brass 

3/4  inch  o.D.  x  1-y  inch 
long  .12-50  microns  pores. 

(d)  Stirrer  I 

Material  : 
Type  : 

Size  : 

Stainless  Steel 

T-line,  electric,  with 
speed  control. 

1/50  H.P.,  20-6000  r.p.m. 
5/16  inch  dia.  shaft. 

(e)  Stirrer  II 

Material  : 
Type  : 

Size  : 

Stainless  Steel 

Lightenin'  Air  Stirrer 
300-3000  r.p.m.  at  30-90  psi 

air  pressure,  5/.l6  inch  dia. 
shaft . 


(f)  Impeller 

Material  : 
Type  : 

Size  : 

Stainless  Steel 

Turbine  type 

1  3/4  inches  dia.,  3  vanes. 

(g)  Reactor  cover 

Material  : 
Size  : 

Stainless  Steel 
9  inch  dia.  x  ^ 

inch  thick 

(h)  Ring  flange 

Material  : 
Size  : 

Stainless  Steel 
6-jj-  inches  i.  d.  , 
-£•  inch  thick. 

9  inches  0. d. , 

Constant  temperature 
bath 

(a)  Vessel 

Material  : 
Size 

Pyrex  glass 

10  inch  dia.  x  12  inches 

(b)  Thermometer 

Material  : 
Range  : 

Mercury-glass 
0-150°  F 

(c)  Heater 

Material  : 
Power  : 

copper 

I5OO  watts 

(CONTINUED) 
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TABLE  B-2 

EQ.UIPMENT  DETAILS  ( CONTINUED) 


(9) 


(10) 


(li) 


(12) 


Equipment 

< 

Constant  temperature  Material  : 

bath  Type 


(d)  Thermoregulator 

Sensitivity: 

(e)  Electronic  Relay 

Material  : 
Type  : 

Capacity  : 

(f)  Stirrer 

Material  : 
Type  : 

Size  : 

Effluent  pump 

Material  : 
Type  : 

Capacity  : 
Plunger  : 

Discharge 
Pressure  : 

Back  pressure 
valve 

Material  : 
Type  : 
Size  : 

Impedance  Cell 

Material  : 
Size  : 

Electrodes: 


Details 


Bimetallic,  American  Inst.  Co. 
+  0.05°  F 


Precision  Scientific  Co. 
on -and- off  control. 

1650  watts 

Stainless  Steel 
Sparkless  'Cenco'  electric 
motor  stirrer 
l/lOO  H.P.,  1000  r.p.m.  , 

40  watts,  -jj-  inch  dia.  shaft, 

1  \  inch  dia.  impeller 

Stainless  Steel 
Positive  displacement,  duplex, 
Milton  Roy  Co. ,  Laboratory 
model 

80  mls/min/side 

-jj-"  dia.  x  1  inch  length, 

96  s/min. 

1000  psia 

Stainless  Steel 
Spring  loaded,  diaphram. 

•5-  inch,  0.6  gpm,  60  psig. 

Lucite 

1  3/4  inches  .length,  ^  inch  i. 
1  \  inches  o.d. 

Platinum,  coated  with  platinum 
black,  26  guage,  \  inch  apart, 
extentending  l/8  inch  in  the 
casing. 


d. 


V  •  ?'  ■ 
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FOR  PRODUCT  WITHDRAWAL 


mT£R  FLOW  RATE  mis -min’ 
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FIGURE  0-2  ^ 


CALIBRATION  OF  ROTAMETER  FOR  WATER  FEED 


ETHYL  ACETATE  FLOW  RATE  mis  *  min" 
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FIGURE  B-3  *  CALIBRATION  OF  ROTAMETER  FOR 


ETHYL  ACETATE 


FEEO 


APPENDIX  C 


EXPERIMENTAL  OBSERVATIONS  AND  METHODS  OF  CALCULATIONS 

C-l  Response  of  the  Reactor  to  a  Step  Function 

The  calibration  of  the  impedance  cell  in  terms  of  the 
concentrations  of  HC.1  solution  versus  its  impedance  as  measured 
in  the  cell  is  given  in  TABLE  (C-l)  and  shown  in  FIGURE  (C-l). 
The  following  relation  was  found  to  hold: 


(C-l) 


IMP  =  Impedance 


ohms 


where 


-1 


CN  =  Concentration  of  HC.1  solution,  gmmoles*ml 

The  calibration  had  good  reproducibility  and  the  impedance 
was  found  to  be  insensitive  to  the  changes  in  flow  rates  encountered 
during  the  experiments.  The  response  of  the  reactor  to  a  sudden 
introduction  of  distilled  water  in  0*.1N  HC.1  solution,  was  followed 
by  the  change  in  the  impedance  of  the  effluent  stream.  The 
experimental  conditions  and  results  for  three  sets,  l-A,  1-B,  l-C, 
are  given  in  TABLES  (C-2)  to  (C-4),  respectively. 


’ 

CELL  IMPEDANCE  ohms 
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FIGURE  C-l  *  CALIBRATION  OF  THE  IMPEDANCE  CELL; 


CONCENTRATION  OF  HCI  SOLUTION  VERSUS  IMPEDANCE 


50.0  S0  3 
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TABLE  C-l 

CALIBRATION  OF  THE  IMPEDANCE  CELL 


Concentration  of  HC1 

solution  CN  S1111110168  x  10^ 
ml 

Impedance 

1.00 

l40 

0.90 

151 

0.80 

.165 

0.70 

185 

0.60 

216 

0.50 

260 

0.40 

330 

0.30 

423 

0.25 

495 

0.20 

608 

0.15 

780 

0.10 

.1.171 

0.06 

1845 

0.03 

3780 

1.00 

l4o 

0.60 

218 

0.25 

493 

0.10 

.1.178 

L 


' 
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TABLE  C-2 


RESPONSE  OF  THE  REACTOR  TO  A  STEP  FUNCTION 


Set  No:  1-A 

Experimental  conditions  during  the  set: 


t 

min 


N  = 

75 

r  .p.m. 

t  = 

45 

min. 

1  = 

44.  5 

mis  •min’”'*' 

¥  = 

97 

gms 

T  = 

80 

°F 

Imp 

CN 

CN 

t  Imp 

CN 

CN 

ohms 

gmmo.lesv-|  n4 

CN(0) 

min  ohms 

gmmoles^i  n4 

CN(0) 

ml 

ml 

0 

l4o 

1.000 

1.000 

24 

221 

0.585 

0.585 

1 

143 

0.966 

0.966 

25 

227 

0.570 

0.570 

2 

144 

0.956 

0.956 

27 

236 

0.549 

0.549 

3 

146 

0.937 

0.937 

28 

240 

0.540 

o.54o 

4 

148 

0.919 

0.919 

29 

245 

0.530 

0.530 

5 

151 

0.895 

0.895 

30 

250 

0. 520 

0.520 

6 

155 

0.865 

0.865 

31 

255 

0. 510 

0.510 

7 

157 

0.851 

0.851 

32 

261 

0.499 

0.499 

8 

160 

0.831 

0.831 

33 

266 

0.490 

0.490 

9 

163 

0.812 

0.812 

34 

271 

0.48.1 

0.48,1 

10 

168 

0.783 

0.783 

36 

282 

0.464 

0. 464 

11 

172 

0.760 

0.760 

37 

289 

0.454 

0.454 

12 

175 

0.745 

0.745 

38 

295 

0.445 

0.445 

13 

178 

0.732 

0.732 

39 

302 

0.436 

0.436 

l4 

181 

0.717 

0.717 

40 

309 

0.426 

0.426 

15 

185 

0.700 

0.700 

42 

322 

o.4io 

0.4.10 

16 

189 

0.684 

0.684 

43 

329 

0.401 

0.401 

17 

192 

0.673 

0.673 

44 

336 

0.3915 

0.3915 

18 

196 

0.659 

0.659 

45 

343 

0.3825 

0.3825 

19 

200 

0.645 

0.645 

46 

350 

o.374o 

0.3740 

20 

205 

0.630 

0.630 

48 

365 

0.3570 

0.3570 

21 

209 

0.617 

0.617 

49 

372 

0.3500 

0.3500 

22 

212 

0.609 

0.609 

51 

389 

0.3320 

0.3320 

23 

217 

0.596 

0.596 

53 

403 

0.3175 

0.3175 

TABLE  (C-2)  continued 


ooo .  r 
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TABLE  C-2  (CONTINUED) 


t 

min 

Imp 

ohms 

CM 

gmolesxl0lt 

ml 

CN 

CN(0) 

55 

422 

0.3010 

O.3OIO 

57 

440 

0.2870 

0.2870 

59 

458 

0.2730 

0.2730 

6i 

477 

0.2620 

0.2620 

63 

497 

0.2490 

0.2490 

65 

520 

0.2370 

0.2370 

67 

543 

0.2250 

0.2250 

69 

568 

0.2150 

O.215O 

71 

591 

0.2060 

0.2060 

73 

6i4 

0.1975 

0.1975 

75 

64o 

O..188O 

O.I88O 

77 

668 

O.I79O 

0.1790 

79 

700 

0.1700 

O.I7OO 

8l 

730 

0.1624 

0.1624 

83 

760 

O..155O 

O.I55O 

85 

793 

0.1480 

0.1480 

87 

827 

0.14.15 

0.14.15 

89 

863 

0.1352 

0.1352 

91 

899 

0.1300 

0.1300 

93 

940 

0.1245 

0.1245 

t 

min 

Imp 

ohms 

CN 

gramoles„in4 

ml 

CN 

cn(o) 

95 

979 

0.1200 

0.1200 

97 

1020 

O.II5O 

0.1150 

99 

IO58 

0.1110 

0.1110 

.101 

1102 

0.1064 

0.1064 

103 

1153 

0.1012 

0.1012 

105 

.1204 

O.O96.I 

0.0961 

109 

1310 

0.0888 

0.0888 

113 

.1428 

O.O8.IO 

0.0810 

117 

1556 

0.0738 

0.0738 

121 

1695 

0.0670 

0.0670 

125 

.1838 

0.0605 

0.0605 

129 

1986 

0.0552 

0.0552 

133 

2.151 

O.O5OI 

0.0501 

137 

2330 

0.0457 

0.0457 

i4i 

2530 

0.04.15 

0.04.15 

145 

2745 

0.0380 

0.0380 

149 

2970 

0.0351 

0.0351 

155 

3350 

0.0320 

0.0320 

158 

3540 

0.0317 

0.0317 

oioz.o 
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RESPONSE  OF 

TABLE  C- 

THE  REACTOR 

-3 

TO  A  STEP  FUNCTION 

Set  No: 

1-B 

Experimental 

conditions  during  the  set: 

N  = 

450 

r  .p.m. 

t  = 

45 

min. 

q.  = 

44. 5 

mis -min' 

-.1 

¥  = 

97 

gms 

T  = 

80 

°F 

t 

Imp 

ON 

CN 

t 

Imp 

CN 

CN 

min 

ohms 

gmmo.les„104  CN(0) 
ml 

min 

ohms 

gmmoles„^A4 

ml 

cn(o) 

0 

i4o 

1.000 

1.000 

55 

4.15 

0.307 

0.307 

1 

140 

1.000 

1.000 

58 

443 

0.284 

0.284 

3 

145 

0.945 

0.945 

6.1 

474 

0.263 

0.263 

5 

.151 

O.895 

0.895 

64 

506 

0.244 

0.244 

7 

.158 

0.843 

0.843 

67 

539 

0.227 

0.227 

9 

165 

0.800 

0.800 

70 

572 

0.2,155 

0.2155 

11 

.171 

0.765 

0.765 

73 

608 

0.200 

0.200 

13 

178. 

5  0.729 

0.729 

76 

647 

0.186 

0.186 

15 

186 

O.697 

0.697 

79 

687 

0.174 

0.174 

17 

.194 

0.666 

0.666 

82 

735 

O.161 

0.161 

19 

202 

0.639 

0.639 

85 

78.1 

O.I5O 

0.150 

2.1 

2.10 

0.615 

0.615 

88 

841 

0.139 

0.139 

23 

2.17 

0.596 

0.596 

91 

895 

0.1305 

0.1305 

25 

225 

0.574 

0.574 

94 

952 

0.123 

0.123 

27 

236 

0.549 

0.549 

97 

.10.10 

0.116 

0.116 

29 

245 

0.530 

0.530 

99 

1055 

0.1113 

0.11,13 

31 

256 

0.508 

0.508 

.104 

1180 

O.O99 

0.099 

33 

266 

0.490 

0.490 

.109 

.1320 

0.088 

0.088 

35 

37 

39 

4.1 

43 

45 

47 

49 

52 


277 

289 

301 

314 

327.5 

34.1 

355 

369 

391 


0.473 

0.454 

0.436 

0.420 

0.402 

0.385 

0.368 

0.3525 

0.329 


0.473 

0.454 

0.436 

0.420 

0.402 

0.385 

0.368 

0.3525 

0.329 


,11 4 
.119 
124 
.129 
134 
139 
144 
149 
154 


i46o 

1635 

.1805 

2008 

2225 

2490 

2780 

3100 

3445 


0.079 

0.0698 

0.0620 

0.0545 

0.0483 

0.0423 

0.0376 

0.0340 

0.0316 


0.079 

0.0698 

0.0620 

0.0545 

0.0483 

0.0423 

0.0376 

0.0340 

0.0316 
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TABLE  C-4 


RESPONSE  OF  THE  REACTOR  TO  A  STEP  FUNCTION 


Set  No:  ,1-C 


Experimental  conditions  during  the  set: 

N  =  55  r.p.m. 

t  =  14.8  min. 

q  =  .135  mls'min 

W  =  97  gms 

T  =  80  °F 


t 

min 

Imp 

ohms 

CN 

gIninoleSrl04 

ml 

CN 

cn(o) 

t 

min 

Imp 

ohms 

CN 

gmmoles„in4 

ml 

CN 

CN(0) 

0 

l40 

1.000 

1.000 

26 

713 

0.167 

O.167 

1 

146 

0.937 

0.937 

27 

760 

0.155 

0.155 

2 

156 

0.858 

0.858 

28 

8.10 

0.145 

0.145 

3 

167 

0.788 

0.788 

29 

863 

0.135 

0.135 

4 

178 

0.732 

0.732 

30 

929 

0.126 

0.126 

5 

.189 

0.685 

0.685 

31 

989 

0.1185 

0.1185 

6 

201 

0.643 

0.643 

32 

.1049 

0.112 

0.112 

7 

214 

0.603 

0.603 

33 

.1.1.10 

0.1055 

0.1055 

8 

227 

0.570 

0.570 

34 

1.180 

0.099 

0.099 

9 

242 

0.536 

0.536 

35 

.1260 

0.0925 

0.0925 

10 

257 

0.506 

0. 506 

36 

1350 

0.086 

0.086 

11 

274 

0.476 

0.476 

37 

.1445 

0.080 

0.080 

12 

293 

0.447 

0.447 

38 

1550 

0.074 

0.074 

13 

314 

0.420 

0.420 

39 

1650 

O.O69 

0.069 

l4 

336 

0.391 

0.391 

40 

.1760 

o.o64 

o.o64 

15 

358 

0.365 

0.365 

4.1 

1870 

0.0595 

0.0595 

16 

380 

0.34.1 

0.34.1 

42 

1990 

0.055 

0.055 

17 

4o4 

0.318 

O.318 

43 

2120 

0.051 

0.051 

18 

430 

0.294 

0.294 

44 

2270 

0.047 

0.047 

19 

458 

0.273 

0.273 

45 

2410 

o.o44 

o.o44 

20 

487 

0.256 

0.256 

46 

2560 

o.o4i 

o.o4i 

21 

518 

0.2375 

0.2375 

47 

2720 

0.0385 

0.0385 

22 

553 

0.221 

0.221 

48 

2900 

0.036 

0.036 

23 

587 

0.212 

0.212 

49 

3100 

0.034 

0.034 

24 

633 

O.I9O 

0.190 

50 

3300 

0.0325 

0.0325 

25 

672 

O.I78 

O.I78 

51 

3500 

0.031 

0.031 

' 

<• 

cl 
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C-2  Sample  Calculations  for  Ester  Conversion 

If,  (a)  Sample  size  of  the  effluent  stream 

(b)  Normality  of  acetic  acid  in  the 
effluent  stream 

(c)  Normality  of  NaOH  used  for  titration 

(d)  mis  of  N/40  NaOH  used  for  titrating 
25  mis  of  effluent  stream 

(e)  mis  of  N/40  NaOH  required  for 
titrating  25  mis  of  water  kept  in 


then , 


contact  with  the  resins 

(f)  mis  of  N/40  NaOH  necessary  to 

titrate  acetic  acid  in  the  25  mis 
of  effluent  stream 


and 


NA  = 


V 


(25) (40) 


=  V  x  10 


-3 


25  mis 

NA 

i/4o 

vs 


0.52 


V 

VS-0.52 


therefore,  Concentration  of  acetic  acid  in 
the  effluent  stream 


Now, 


!e  ?e/: 


m« 


^e  +  %r 
leSe 


q  Me 


A, 


o 


V  x  10"6  gmmoles •ml--'- 


' 
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where, 

Cj_  =  Influent  concentration  of  ethyl  acetate  gmmoles*ml-1 
M0  =  Molecular  weight  of  ethyl  acetate 


=  88.11 

qg  =  Flow  rate  of  ethyl  acetates 
q^  =  Flow  rate  of  water 
q  =  Total  influent  flow  rate 

=  S*  +  S* 

=  Density  of  ethyl  acetate 
=  0.892 


mis* min-1 
mls*min-1 
mis •min-1 

gms*ml_1 


Since  one  mole  of  ethyl  acetate  hydrolysed  yields  one 
mole  of  acetic  acid, 

°lo  conversion  of  ethyl  acetate  =  X(.100) 

A 

=  —  x  100 

Ci 


detailed  calculations  for  Set  3-A 


Run  No.  3-A-4 


ci 


9-eSe 

qMe 


1 

(5) 

(0.892) 

(100) 

(88.11) 

=  5.0732  x  10-4 

=  6.4i 


gmmoles*ml 


-1 


Y 


mis 
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An  —  V  X  10 


-6 


=  6.4l  x  10 


-6 


and  X(.100)  =  (,6-lfl  x  10 — x  .100 

(5.0732  x  .lCT4) 


gmmo.les*ml 


-1 


=  1.26 

£i  =  5.0732  x  10~4 

^•o  6 . 4.1  x  10“^ 

=  79.15 

Similar  results  for  the  other  runs  in  set  3-A  are  given  in 
TABLE  (C-8).  Also,  the  above  calculations  for  the  experimental 
runs  made  to  evaluate  the  significance  of  external  film  diffusion 
are  shown  in  TABLES  (C-5)  to  (C-7). 

C  • 

The  functional  relation  between  _J-  and  q  is  given  by 

-A-o 

(equation  (IV.18)  ): 


1  ♦  3 


Q 


This  is  an  equation  of  a  straight  line  with  slope  equal 

to  i  and  intercept  of  unity  with  the  vertical  axis.  The  best 
Q  C- 

fit  of  the  experimental  data  of  _i  vs  q  was  taken  as  the  line 

■A-0 

which  passes  through  (0,.l)  and  gives  the  least  sum  of  squares  of 
the  deviations  of  the  individual  observations  from  the  line.  It 
can  be  shown  that  for 


( 


' 


TABLE  C-5 


SIGNIFICANCE  OF  EXTERNAL  FILM  DIFFUSION 


Set  No:  2-A 

Experimental  conditions  during  the  set: 


T  = 

80 

0  y 

W  = 

95 

gms 

R  = 

0.00283 

cms 

38.4 

mls*min-^ 

^e  = 

1.6 

mls*min“l 

q.  = 

4o.o 

mls*min~l 

ci  = 

5.0732  x  io_if 

gmmoles.ml'^ 

t  = 

50 

min 

Run  No. 

N 

V* 

Ad 

X(lOO) 

Ci 

r  .p.m. 

mis 

gramoles,rlrt6 

ml 

°/o  conversion 

Ao 

A-l 

50 

3.91 

3.91 

0.77 

129.75 

A- 2 

100 

14.01 

14.01 

2.76 

36.21 

A- 3 

125 

32.12 

32.12 

6.33 

15.79 

A- 4 

150 

50.22 

50.22 

9.90 

10.10 

A- 5 

200 

52.46 

52.46 

10.34 

9.67 

A- 6 

275 

53.77 

53.77 

10.60 

9.43 

A- 7 

325 

52.93 

52.93 

10.44 

9.58 

A- 8 

4oo 

53.64 

53.64 

10.57 

9.46 

A- 9 

475 

53.48 

53-48 

10.54 

9.49 

*  V  =  mis  of  N/40  NaOH  necessary  to  titrate  acetic  acid  in  25  mis 
of  the  effluent  stream. 
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TABLE  C-6 


SIGNIFICANCE  OF  EXTERNAL  FILM  DIFFUSION 


Set  No.  2-B 


Experimental  conditions  during  the  set: 


T  = 

80 

0  y 

W  = 

48 

gms 

R  = 

0.00283 

cms 

q.  = 

w 

38.4 

mis*  min"1- 

q.  = 

e 

1.6 

mis -min"1 

q.  = 

4o.o 

mis -min"1 

ci  = 

5.0732  x  10_1+ 

gmmoles’inl”-*- 

t  = 

50 

min 

Run  No. 

N 

V 

^0 

x(ioo) 

Ci 

r  .p.m. 

mis 

gmmolesvn06  conversion 

Ao 

ml 

B-l 

75 

2. 

58 

2.58 

0.51 

196.64 

B-2 

125 

6 . 

88 

6.88 

1.36 

73.74 

B-3 

150 

16. 

86 

16 . 86 

3.32 

30.09 

b-4 

200 

28. 

47 

28.47 

5.61 

17.82 

B-5 

320 

27. 

92 

27.92 

5.50 

18.17 

b-6 

420 

28. 

6.1 

28.61 

5.64 

17.33 

B-7 

500 

28. 

53 

28.53 

5.62 

17.78 
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TABLE  C-7 

SIGNIFICANCE  OF  EXTERNAL  FILM  DIFFUSION 

Set  No.  2-C 

Experimental  conditions  during  the  set: 


T 

=  100 

0  F 

W 

= 

48 

gms 

R 

= 

0.00283 

cms 

q. 

w 

= 

38.4 

mls*min"^ 

q. 

e 

= 

1.6 

mis -min-1 

q. 

- 

4o.o 

mis  •min”"1' 

ci 

= 

5.0732  x  10* ^ 

gmmoles'ml"^ 

t 

= 

50 

min 

Run  No. 

N 

V 

Ad 

X(.100) 

Ci 

r  .p.m. 

mis 

gmmoles,r^n6 

ml 

i  conversion 

Ac 

C-.1 

100 

17-37 

17.37 

3-42 

29 . 2.1 

C-2 

200 

68. 46 

68.46 

13.50 

7.4.1 

C-3 

350 

67.46 

67.46 

13.30 

7.52 

c-4 

500 

68.28 

68.28 

13.46 

7.43 

✓ 
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n  2 


X 


n  =  number  of  runs  in  a  set 

to  be  minimum,  the  slope  of  the  line  is 


Thus  for  set  3- A  we  have  the  following  values  (See  TABLE  (C-8)) 


Run  No. 

ci 

q. 

c  • 

(-)l 

Ao 

(caY 

k 

Ao 

\Ac J 

3-A-l 

31.43 

4o 

1257.2 

l600 

3- A- 2 

47.10 

60 

2826 . 0 

3600 

3- A- 3 

62.40 

80 

4992.0 

64oo 

3-A-4 

79.15 

100 

7915.0 

10000 

n  =  4 

220.08 

280 

16990.2 

21600 

• 

1 

(16990.2)  - 

(280) 

•  • 

Q 

(21600) 

=  0.7736 


Q 


1.2926 
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For  each  set,  the  values  of  Q  are  obtained  in  this 
manner.  Thus  at  80°  F  we  have  the  following  values  (See  TABLES 
(C-8)  to  (C-ll) : 


Set  No. 

W 

R 

Q 

gms 

cms 

mis -min" 

3-A 

28 

0.03045 

1.2926 

3-B 

42 

0.03045 

1.9466 

3-C 

51 

0.01268 

2.5289 

3-D 

68 

O.OO283 

3.4428 

The  equation 


q  -  (wR-coth(wR)  -  l) 

R2^ 


wher  e , 


k 

D 


(IV. 15) 


was  then  used  to  evaluate  the  parameters  k  and  D  from  the 
experimental  data  shown  above  -  both,  k  and  D,  being  independent 
of  W  and  R.  The  equation  has  two  unknowns,  w  and  D,  and 
there  are  four  sets  of  observations,  Hence,  the  'best'  values  of 
w  and  D  were  obtained  by  the  following  trial-and-error- 
calculations  procedure. 

(1)  Assume  a  reasonable  but  arbitary  value  of  w. 

(2)  Using  equation  (IV.  .15)  calculate  D  for  each 
set  of  experimental  values  of  Q,  W  and  R. 

This  gives,  say,  D^,  Dg,  and  D^. 


‘ 

. 
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(3)  Calculate  the  mean  D  : 

4 

I  Di 

1 


(4)  Calculate  the  standard  deviation  of  D^’s  with 
respect  to  E^,  using  the  following  formula 

Vk  2 

-  dj 

1 

s.D.  = - 

4-1 

(5)  Change  the  value  of  w  and  repeat  steps  (2)  to  (5). 

(6)  Select  the  range  of  w  within  which  the  minimum 
standard  deviation  as  calculated  in  step  (4)  lies. 

(7)  Now  vary  w  hy  much  smaller  values  and  repeat  steps 
(l)  to  (6). 

Thus  progressively  the  range  of  w  was  narrowed  down 
and  smaller  jumps  were  given  in  varying  w  .  In  the  final  set  of 
calculations  the  successive  values  of  w  differed  hy  only  0.01. 
Eventually,  the  value  of  w  which  gave  the  minimum  standard 
deviation  of  D^'s  with  respect  to  the  mean,  E^  (Step-4),  was 
chosen  as  the  best  value  and  the  corresponding  E^  as  the  best 
value  for  diffusivity,  D. 

Using  the  'best'  values  of  w  and  D,  Q,'s  were  calculated 
using  equation  (lV-15)  and  were  denoted  QC.  The  final  values  for 
80°  F  were  as  follows: 
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Set  No. 

Q 

D  x  105 

QC 

mis •min-1 

cm2,  min"1 

mis* min 

3-A 

1.2926 

2.9805 

1 . 3021 

3-B 

1 . 9466 

2.9719 

1.9532 

3-C 

2.5289 

2.9978 

2.5248 

3-D 

3.4428 

3.0366 

3*4.119 

mean 

2.9967 

w  = 

36.52 

D  - 

2.9967  x  10" 5 

2  •  -1 

cm  *min  ^ 

k  = 

w2D 

= 

3.997  x  10“ 2 

min"1 

The  Relative  Standard 

Deviation  (RSD)  of 

the  above 

was  calculated  by  using  the  following  expression 

4  2 

1 

2 

7  (Di  -  ig 

RSD  = 

.1 

X 

100 

_  (4  -  1)%2 

= 

0.67  i 

Results  of  similar  calculations  for  100° 

F  &  120°  F 

-1 


are  shown  in  TABLE  (V-2) . 
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TABLE  C-8 

EVALUATION  OF  k  AND  D 


Set  No,  3-A 

Experimental  conditions  during  the  set: 


T  = 

80 

o  y 

W  - 

28 

gms 

R  = 

0,03045 

cms 

N  = 

600 

r  ,p  .m. 

Gi  = 

5.0732  x  .10' 

gmmoles  •m.l-^ 

Run  No. 

^e 

mls/min 

% 

mls/min 

q  t 

mls/min  min. 

V  Ao 

mis  gmmoles  ^2_n6 
ml 

Gi 

a0  i 

X(.100) 

conversion 

A-l 

2.0 

38 

4o  50 

16.. 1 4  16.  l4 

31.43 

3.18 

A- 2 

3.0 

57 

60  33.33 

10.77  10.77 

47.10 

2.12 

A- 3’ 

4.0 

76 

80  25 

8.13  8.13 

62.40 

1.60 

A- 4 

5-0 

95 

100  20 

6 . 4.1  6 . 4l 

79.15 

1.26 

The  least  square  fit  of  the  equation  (lV.l8)  gives 

-  =  0.7736 

Q 

Q  =  I.2926  ml  s  •  min 
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TABLE  C-9 

EVALUATION  OF  k  AND  D 


Set  No.  3-B 

Experimental  conditions  during  the  set: 


T  = 

80 

°  F 

W  = 

42 

gms 

R  = 

0.03045 

cms 

N  = 

600 

r  .p.m. 

Ci  = 

5.0732  x  10"4 

gmmoles*ml" 

1 

Run  No. 

A 

mls/min 

X 

mls/min 

1 

mls/min 

t  V 

min.  mis 

Ao 

a5Solesxl06 

ml 

cx 

Ao  i 

X(.100) 

conversion 

B-l 

1.5 

28.5 

30 

66.67  31.53 

31.53 

16 . 09 

6.21 

B-2 

2.5 

47.5 

50 

4o.oo  20.08 

20.08 

25.26 

3.96 

B-3 

3.5 

66. 5 

70 

28.57  14.42 

14.42 

35.18 

2.84 

B-4 

4.5 

85.5 

90 

22.22  10.67 

10.67 

^7.55 

2.10 

B-5 

5.5 

104.5 

110 

18.17  8.57 

8.57 

59.20 

I.69 

The  least  square  fit  of  the  equation  (lV..l8)  gives 

|  =  °'5137 

.*.  Q  =  1.9466  mls’min 


< 


> 
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TABLE  C-10 

EVALUATION  OF  k  AND  D 


Set  No.  3-C 

Experimental  conditions  during  the  set: 


T  = 

80 

0  -p 

W  = 

51 

gms 

R  = 

0.01268 

cms 

N  - 

6  00 

r  .p.m. 

Ci  - 

5.0732  x  10" 

gmnioles/ml 

Run  No. 

mls/min 

q. 

w 

mls/min 

q  t 

mls/min  min. 

V  Ao 

mis  gmmoles^^ 
ml 

Ci 

A0  * 

x(.100) 

conversion 

C-l 

1.75 

33.25 

35  57.14 

32.71  32.71 

I5.5I 

6.45 

C-2 

3.0 

57 

60  33.33 

20.25  20.25 

25.05 

3.99 

C-3 

4.0 

76 

80  25 

15.71  15.71 

32.29 

3.10 

c-4 

5-0 

95 

100  20 

12.56  12.56 

40.39 

2.4.1 

The  least  square 

fit  of  e 

quation  (IV..18) 

gives 

1 

Q 

0.3954 

• 

•  • 

Q  = 

2.5289 

,  .  -1 
mis*  mm 
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TABLE  C-ll 

EVALUATION  OF  k  AND  D 


Set  No.  3-B 

Experimental  conditions  during  the  set: 


T  = 

80 

0  F 

w  = 

68 

gms 

R  = 

0.00283 

cms 

N  = 

600 

r  ,p.m. 

Ci  = 

5.0732  x  10_1+ 

gmmoles- 

■ml”'1 

Run  No. 

^e 

V 

q  tV 

A 

0 

Ci 

x(ioo) 

mls/min 

mls/min 

mls/min  min.  mis  gmmoles^.^6  ac  °1° 

ml 

conversi 

D-l 

1.0 

19 

20  100  75.59 

75.59 

6.71 

14.90 

D-2 

2.0 

38 

4o  50  38.91 

38.91 

13.04 

7.67 

D-3 

3.25 

6.1.75 

65  30.77  25.99 

25.99 

19.52 

5.12 

d-4 

4.5 

85.5 

90  22.22  18.94 

18.94 

26.79 

3.73 

D-5 

6 . 0 

114 

120  16.67  14.02 

l4. 02 

36.19 

2.76 

The  least  square 

fit  of  equation  (IV.18)  gives 

1 

0, 


0.2905 


Q 


3.4428 


mis*  min 


08 
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TABLE  C-12 

EVALUATION  OF  k  AND  D 


Set  No.  4-A 


Experimental  conditions  during  the  set: 

T  =  100 

W  =  80 

R  =  0.03045 

N  =  600 

C±  =  4.0586  x  10'4 


o  p 

gms 

cms 


r  .p.m. 
gmmoles*ml 


q. 

1 

q. 

t 

V 

A 

n 

x(ioo) 

Run  No. 

e 

w 

c. 

mls/min 

mls/min 

mls/min 

min. 

mis  gmmoles„106 
ml 

1 

Ad 

$  conversion 

A-l 

1.2 

28.8 

30 

66.67 

97.23 

97.23 

4.17 

23.98 

A- 2 

2.4 

57.6 

60 

33.33 

56.56 

56.56 

7.18 

13.93 

A- 3 

3.2 

76.8 

80 

25 

42.29 

42.29 

9.60 

10.42 

A- 4 

4.0 

96.0 

100 

20 

34.09 

34.09 

II.91 

8.40 

The  least  square  fit  of  equation  (IV.18) 


gives 


i  =  O.IO74 

Q 

Q  =  9.3.IOO  mis*  min 


' 


. 
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TABLE  C-13 

EVALUATION  OF  k  AND  D 

Set  No.  4-B 


Experimental 

conditions 

during  the  set: 

T  = 

100 

0  F 

W  = 

43 

gms 

R  = 

0.01268 

cms 

N  = 

600 

r  .p.m. 

C.  = 

1 

4.0586  x  10'4 

gmmoles*ml~ 

1 

1 

Run  No.  e 

mls/min 

q. 

w 

mls/min 

q  tV 

mls/min  min.  mis 

A 

0 

gmmolesvln6 

ml 

C. 

r  % 

Ao 

x(ioo) 

conversion 

B-l  1.2 

28.8 

30  66.67  59.70 

59.70 

6.80 

14.71 

B-2  2.0 

48.0 

50  40.00  39.93 

39.93 

10.16 

9.84 

B-3  2.8 

67.2 

70  28.57  29.93 

29.93 

13.56 

7.37 

B-4  3.6 

86.4 

90  22.22  23.52 

23.52 

17.26 

5.79 

The  least  square  fit  of  equation  (lV.l8)  gives 


i  =  0.l8l4 

Q 


Q  = 


mis*  min 


-1 


•  « 


5.5137 
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TABLE  C-l4 

EVALUATION  OF  k  AND  D 


Set  No.  4-C 


Experimental  conditions 

during  the  set: 

T  = 

100 

0  F 

W  = 

70 

gms 

R  = 

O.OO283 

cms 

N  = 

600 

r  .p.m. 

C.  = 

l 

4.0586  x  10" 

•4 

gmmoles*ml‘ 

-1 

Run  No. 

q  t 

V 

A 

0 

x(ioo) 

Ci 

mls/min 

mls/min 

mls/min  min. 

mis  gmmoles„1Pl6 
ml 

°Jo  convers>ion 

C-l  1.6 

38.4 

4o  50 

75.90 

75.90 

5.35  18.69 

C-2  2.4 

57.6 

60  33.33 

51.44 

51.44 

7.89  12.67 

C-3  3.2 

76.8 

80  25 

4o.66 

4o.66 

9.98  10.02 

c-4  4.0 

96.0 

.100  20 

32.81 

32.81  12.37  8.08 

The  least  square 

fit  of  e 

quation  (lV.l8) 

gives 

1 

Q 

O.II3I 

m 

•  • 

Q  = 

8.84,17 

mls.min--1- 
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TABLE  C-15 

EVALUATION  OF  k  AND  D 


Set  No.  5-A 


Experimental 

conditions 

during  the  set: 

T  = 

120 

0  F 

W  = 

80 

gms 

R  = 

0.03045 

cms 

N  = 

600 

r  .p.m. 

I  = 

4.0586  x  .10' 

-4 

gmmoles*ml' 

-1 

qe 

% 

q  t 

V 

Ao 

x(.100) 

Run  No. 

^i 

mls/min 

mls/min 

mls/min  min. 

mis  gmmoles„-j^6 
ml 

A0  %  conversion 

A-l  1.2 

28.8 

30  66.67  I7O.O9 

170.09 

2.39  4,1.84 

A- 2  2.0 

48.0 

50  4o 

128.05 

128.05 

3.17  31.54 

A- 3  3.2 

76.8 

80  25 

86.42 

86.42 

4.70  21.28 

A- 4  4.4 

105.6 

110  18.17 

64.38 

64.38 

6.30  15.87 

The  least  square  fit  of  equation  (IV.18)  gives 

i  =  0.0470 

0, 

Q  21.2787  mis •min 


, 

f 
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TABLE  C-l6 

EVALUATION  OF  k  AND  D 


Set  No.  5-B 

Experimental  conditions  during  the  set: 


T  = 

120 

0  y 

W  = 

43 

gms 

R  = 

0.01268 

cms 

N  = 

600 

r  .p.m. 

ci  = 

4.0586  x  10-^ 

gmmoles*ml‘ 

-1 

Run  No. 

% 

1 

t 

V 

Ad 

Ci 

x(.100) 

mls/min 

mls/min 

mls/min 

min. 

mis  gmmoles^-i  n6 

Ao  * 

convers 

ml 

B- 

•1 

1.2 

28.8 

30 

66.67 

124.24 

124.24 

3.27 

30.58 

B- 

•2 

2.0 

48.0 

50 

40 

83.29 

83.29 

4.87 

20.53 

B- 

■3 

2.8 

67.2 

70 

28.57 

6.1.38 

6.1.38 

6.6.1 

15.13 

B- 

■4 

3.6 

86.4 

90 

22.22 

50.43 

50.43 

8.05 

12.42 

The  least  square  fit  of  equation  (IV.18)  gives 


^  =  O.O786 


Q  = 


mis*  min 


-1 


12.7251 
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TABLE  C-17 

EVALUATION  OF  k  AND  D 


Set  No.  5-C 

Experimental  conditions  during  the  set: 

T  =  120 

W  = 

R  = 

N  = 


Ci  = 


70 

0.00283 

600 

4.0586  x  10“' ^ 


0  F 
gms 
cms 

r  .p.m. 

gmmoles*ml 


-1 


Run  No. 

^e 

mls/min 

q. 

w 

mls/min 

q. 

mls/min 

t 

min. 

V  Ao 

mis  gmmoles,r^n6 
ml 

X(100) 

Li 

_4o  $  c  on' ver  s  i  on 

C-l 

1.6 

38.4 

4o 

50 

142.23 

142.23 

2.85  25.09 

C-2 

2.4 

57.6 

60 

33.33 

103.79 

103.79 

3.91  25.57 

0-3 

3.2 

76.8 

80 

25 

80.53 

80.53 

5.04  19.84 

c-4 

4.0 

96.0 

100 

20 

70.59 

70.59 

5.75  17.39 

The  least  square 

fit  of  e 

quation  (lV.l8) 

gives: 

1 

Q 

0.0485 

• 

•  • 

Q  = 

20.6343 

mis  ♦min"  -'- 

' 

